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CHAPTER  I 


INTRODUCTION 

This  work  considers  the  development  of  new  materials  for  applica¬ 
tions  in  infrared  nonlinear  optics.  The  investigated  crystals  belong  to 
the  Il-IV-Vg  class  of  ternary  semiconductors  with  chalcopyrite  struc¬ 
ture.  We  discuss  growth  technique,  materials  quality  problems,  ab¬ 
sorption  mechanisms,  semiconducting  properties,  and  linear  and  non¬ 
linear  optical  properties  for  two  of  the  compounds:  CdGeASg  and 
CdGePp  . 

In  general  nonlinear  crystals  are  used  to  generate  coherent 
radiation  at  wavelengths  where  lasers  are  not  available.  The  process 
can  be  second  harmonic  generation  ( SHG ) ,  mixing,  or  parametric 
oscillation. 

Since  nonlinear  optics  require  high  optical  power  densities 
the  material  quality  problems  are  extremely  important.  Alto,  the 
crystals  must  be  phasematchable  and  have  a  high  nonlinear  figure  of 
merit  to  give  good  nonlinear  interaction  efficiencies.  Several  use¬ 
ful  crystals  exist  for  the  visible  spectral  region.  In  the  infrared, 
however,  there  is  a  lack  of  high  quality  nonlinear  crystals. 

As  a  result  of  this  study,  CdGeASg  has  been  shown  to  have  the 
highest  figure  of  merit  of  any  presently  available  infrared  nonlinear 
material.  With  a  positive  birefringence  of  0.1,  CdGeASg  is  phase¬ 
matchable  for  SHG  and  parametric  interactions  over  most  of  its  trans- 
pnency  range  from  2.U  to  18  pm.  By  SHG  of  a  TEA  COg  laser  we  have 
observed  a  doubling  efficiency  of  2  percent. 
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To  improve  the  understanding  of  the  optical  absorption  mechanisms 
we  have  measured  the  semiconductor  properties  and  determined  the  band 
structure  near  the  Brillouin  zone  center.  We  have  identified  one 
of  the  absorption  mechanisms  in  CdGeAs2  as  due  to  intraband  transitions 
between  the  split  valence  bands.  This  absorption  is  avoided  by 
reducing  the  hole  concentration  through  compensation  or  doping. 

The  II-IV-V2  chalcopyrite  compounds  were  first  synthesized  in 
the  late  nineteen  fifties.1'2  Over  the  past  ten  years  Russian  workers 
have  reported  extensive  work  on  the  compounds.'"’  We  had  therefore 
information  about  melting  points  and  several  attempted  growth  methods 
when  the  chalcopyrite  growth  program  started  at  Stanford  in  the  fall 
of  1969.  The  important  properties  for  nonlinear  optics  such  as  dis¬ 
persion,  birefringence ,  and  nonlinear  coefficients  were  not  known, 
but  from  the  close  resemblance  with  the  III-V  compounds  we  assumed  that 
the  nonlinear  coefficients  were  large.  From  the  crystal  structure  it 
was  clear  that  the  crystals  would  be  birefringent  because  of  the 
tetragonal  distortion  of  the  unit  cell.  We  did  not  know,  however, 
if  the  birefringence  would  be  sufficient  to  satisfy  the  phasematching 
condition.  Today,  due  to  the  recent  work  at  Stanford  and  the  Bell 
Telephone  Laboratories,  the  potential  of  the  Il-IV-Vg  chalcopyrite 
compounds  for  nonlinear  optics  is  fully  confirmed. 

The  long-term  goal  of  this  research  is  tc  develop  chalcopyrite 
crystals  of  high  enough  optical  quality  for  the  construction  of  a 
parametric  oscillator  tunable  over  the  infrared  frequency  range.  In 
Chapter  VII  we  discuss  the  potential  use  of  tunable  infrared  radiation 


2 


to  detect  air  pollutants  and  compare  the  sensitivity  with  alternative 
laser  detection  schemes. 


t 


1 


CHAPTER  II 


SELECTION  OF  CRYSTALS  FOR  NONLINEAR  OPTICS 


A.  GENERAL  CONSIDERATION 


This  chapter  reviews  briefly  the  crystal  parameters  ol  .•  .portance 
for  nonlinear  optics  and  introduces  the  properties  of  the  chalcopyrite 
compounds.  To  evaluate  nonlinear  crystu.  s  it  is  necessary  to  considei 
the  following  properties: 


1.  Large  nonlinear  coefficient. 

2.  Sufficient  birefringence  for  phasematching. 

3.  Crystal  transparency  range, 
k.  Optical  quality. 


The  second  order  nonlinear  coefficient  d^  is  defined  by 


,+#5 


Pil’(tU5)  =  dijkEj(uxL)Ek(a52} 


(2.1) 


with  summation  over  repeated  indices.  In  this  equation  P^'u^)  is  the 
nonlinear  polarization  generated  by  the  electric  fields  EjCu^)  and 
Ek(u^).  The  symmetry  of  the  d^  tensor  follows  from  the  crystal 
point  group  symmetry.^  For  the  crystal  to  have  a  nonvanishing  nonlinear 
coefficient  it  cannot  have  a  center  of  symmetry.  In  materials  where  the 
nonlinear  coefficient  has  not  been  measured ,  the  magnitude  can  be 
estimated  by  Miller's  rule^  if  e  indices  of  refraction  are  known.  The 
Miller's  rule  relates  the  indices  of  efraction  n  to  the  nonlinear 
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coefficient  by 


dijk  -  -  ‘H-j  '  D(^  -  1)  ,  (2-2) 

where  6...  is  Lhe  Miller’s  constant.  For  most  materials,  6.,,  is 
1JK  ’  ljk 

constant  to  within  a  factor  of  two  while  the  nonlinear  coefficient  can 
vary  over  almost  four  orders  of  magnitude.  Miller's  rule  works  very 
well  especially  for  crystals  with  the  same  point  group  symmetry. 

The  nonlinear  polarization  drives  an  electric  field  at  .  In 

order  for  the  driving  polarization  to  stay  in  phase  with  the  resulting 
electric  field  over  the  whole  crystal  length  the  phasematchin0  condition 


k5  =  kl  +  k2 


(2.3) 


mus 


t  be  satisfied.  Since  k  =  nui/c  ,  we  can  rewrite  the  equation  as 


n^  =  njCOj  +  ngOJg 


(2.4) 


In  addition,  energy  conservation  requires 


*3 


(2.5) 


In  crystals  we  accomplish  phasematching  by  using  the  birefringence 
to  compensate  for  the  dispersion  in  the  indices  of  refraction.  The 
condition  that  the  crystals  should  have  sufficient  birefringence  to  achieve 
phasematching  sets  a  very  important  limit  on  the  number  of  crystals 
useful  for  nonlinear  optics.  Degenerate  parametric  oscillators  or  second 
harmonic  generation  (SHG)  require  the  largest  birefringence.  The 


phasemacching  angle,  for  SHG  is  approximately  given  by  (cf.  Appendix  A ) 

sin  ft)  VD/|B  |  (2.6) 

for  type  I  phasematching  and 

sin  1:53  V^D/I B  (2.7) 

for  type  II  phasematching.  Here  B  -  ng  -  n^  is  the  birefringence 
and  D  =  n^“  -  n^1  is  the  dispersion.  The  crystals  must  therefore  have 
a  birefringence  larger  than  the  dispersion  between  the  second  harmonic 
and  the  fundamental  frequency. 

The  transmission  range  of  crystals  is  determined  by  the  bandgap 
absorption  at  short  wavelengths  and  by  two-phonon  absorption  at  long 
wavelengths.  Since  efficient  nonlinear  interactions  require  high 
optical  power  densities,  the  absorption  within  the  transparency  range 
must  be  at  an  absolute  minimum  to  avoid  crystal  damage  or  breaking  of 
the  phasematching  condition.  For  reference,  in  a  high  optical  quality 

-1  9 

infrared  material  such  as  GaAs  the  absorption  is  0.02  cm  ,  and  the 

2  10 

optical  damage  threshold  is  approximately  100  MW/cm  .  Nonlinear 
absorption  is  also  important  and  may  even  limit  the  maximum  efficiency 
of  the  nonlinear  interaction.  We  will  discuss  this  in  Chapter  VI. 

B.  ESTIMATION  OF  THE  OPTICAL  PROPERTIES  OF  THE  Il-lV-Vg  CHALC0PYRITE 
COMPOUNDS 

The  III-V  semiconductors  have  large  nonlinear  coefficients,  but 
due  to  the  cubic  symmetry  they  have  no  birefringence  and  phasematching 
is  therefore  not  possible.  If  we  replace  the  cation  in  the  III-V 
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compounds  with  atoms  from  the  second  and  fourth  column  of  the  periodic 
table,  we  generate  a  new  class  of  crystals  known  as  the  II-IV-V^ 
chalcopyrite  compounds.  These  have  very  similar  optical  properties 
but  with  the  important  addition  that  they  are  birefringent .  The 

chalcopyrite  crystal  structure  is  tetragonal  and  the  point  group 

r  4 

symmetry  of  42m  allows  both  type  1  and  type  II  phasematching. 

Figure  1  shows  the  atomic  arrangement  in  the  chalcopyrite  unit 

12 

cell.  The  lattice  parameters  are  a  and  c  .  Each  cell  contains 
four  II-IV-V^  formulas.  If  the  II  and  IV  atoms  were  indistinguishable, 
the  chalcopyrite  structure  would  reduce  to  a  compressed  sphalerite 
structure.  However,  because  of  the  ordered  arrangement  of  the  II  and 
IV  atoms,  the  chalcopyrite  unit  cell  is  twice  the  sphalerite  cell.  The 
II  and  IV  cations  have  different  covalent  radii,  and  the  anion  V  is 
therefore  not  exactly  in  the  centtr  cf  the  cation  tetrahedra.  It  is 
positioned  a  distance  a(x  -  1/4)  off  the  center  closer  to  the  cations 
with  the  smallest  radii. ^  We  can  characterize  the  chalcopyrite 
structure  by  two  parameters  r  and  a  where  t  =  2  -  c/a  is  the 
tetragonal  compression  and  a  =  4x  -  1  gives  the  position  of  the 
anion. 

The  growth  of  high  optical  quality  material  presents  several 
problems.  For  maximum  progress  we  decided  to  concentrate  on  a  few 
compounds.  In  choosing  which  compounds  to  grow,  the  important  con¬ 
siderations  were  transmission  range,  birefringence,  and  magnitude  of 
the  nonlinear  coefficients.  When  we  started  the  growth,  the  only 
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useful  information  available  was  the  crystal  structure  and  the  bandgap 
frequency.  The  other  important  properties  such  as  the  infrared  cut¬ 
off  frequency,  the  nonlinear  coefficients,  and  the  birefringence  had 
not  been  measured  e  therefore  used  the  close  resemblance  with  the 
III-V  analogs  to  predict  the  properties.  For  reference  the  III-V 
analogs  are  1  sted  in  Table  1.  As  pointed  out,  two-phonon  absorption 
determines  tht  infrared  cut-off  frequency,  and  the  frequency  is 
approximately  twice  the  transverse  optical  mode  frequency  at  the 
Brillouin  zone  center.  When  we  assume  the  same  force  constants 
for  the  chalcopyr ite  bonds  as  for  the  III-V  analogs,  we  estimate 
the  transverse  optical  mode  frequency  to  be  equal  to  the  largest 
transverse  optical  mode  frequency  of  the  III-V  analogs.  This  leads 
to  predicted  cut-off  wavelengths  that  agree  with  recent  measurements 
to  better  than  one  micron.  For  example,  for  CdGeAsg  we  predicted 
a  cut-off  at  18.5  pm  and  measured  18  pm  as  shown  in  Fig.  2.  The 
strength  of  the  two  phonon  absorption  in  CdGeAs2  is  about  20  cm”1 
at  room  temperature.  There  is  also  some  weak  chree  phonon  absorption 
of  approximately  0.3  cm”1. 

The  chalcopyrites  have  only  one  independent  nonlinear  coefficient 
since  =  according  to  the  Kleinman  symmetry  condition.1^ 

By  applying  the  Miller's  rule  in  Eq.  (2.2),  we  estimate  the  magnitude 
using  the  measured  Miller's  constants  for  the  III-V  analogs.  The 
agreement  with  the  recent  experimental  values  is  within  a  factor  of 
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OPTICAL  PROPERTIES  0?  THE  CKALCOPYRITES 


Two  phonon 
absorption 


FIG,  2--Tvo  and  three  phonon  absorption  in  CdGeAs 


two.  For  compounds  with  unknown  average  index  of  refraction  we  use 
Eq .  (2.10)  to  estimate  the  index. 

The  property  most  difficult  to  estimate  is  the  birefringence.  The 
c-axis  compression  of  the  III-V  compounds  results  in  a  negative 
birefringence,  so  it  was  surprising  that  mo3t  of  '  t«e  chalcopyrites 
turned  out  to  be  positive  birefringent .  Initially,  vc  assumed  that 
the  crystals  rfith  the  largest  tetragonal  distortion  had  also  the  largest 
birefringence.  We  have  shown  later,  however,  that  there  is  no  such 
correlation.  For  example,  both  CdGeAsg  and  CdGePg  have  large  tetragonal 
distortion,  but  only  CdGeAsg  has  large  enough  birefringence  for  phase¬ 
matching  . 

It  is  interesting  to  derive  expressions  for  the  indices  of  refraction 
and  the  birefringence  for  :he  chalcopyrites  by  assigning  bond  polarizabilities 
to  the  II-V  and  the  III-V  bonds.  With  a  linear  polarizability  tensor 
a  ,  we  write  the  induced  polarization  of  a  bond  as 

p(t)  =  e0a-E(t)  •  (2.8) 


We  assume  uniaxial  bond  symmetry  and  write  the  transverse  and  the 
longitudinal  linear  bond  polarizabilities  as  and  ct jj  .  We  can 

then  calculate  the  susceptibility  tensor  by  summing  over  the  thirty-cwo 
bonds  within  the  chalcopyrite  unit  cell.  Defining  an  average  bond 
polarizability  as  a  =  j(2aj^  +  Ct||)  3nd  the  anisotropy  as  7  =  (a|j  -  a^). 
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we  follow  Ch.rula1^  and  obtain  for  the  susceptibility  tensor  that 


X  =  €  -  1  = 


l .6  16  f 

7  <“aC  +  °B0>i+  ~  pAC(T  +  + 


/l  0  0  \ 

cr)l(o  1  0  ), 
J\o  0  -2/ 

(2.9) 


where  AC  and  BC  refer  to  the  bonds  to  the  largest  and  smallest 
cation  and  V  is  the  unit  cell  volume.  We  assume  the  average  bond 
polarizability  is  approximately  the  same  as  for  the  IIJ-V  analogs  and 
calculate  the  average  susceptibility  for  the  chalcopyrites  using  the 


equation 


V  vk 

-  x  +  —  \ 

V  3  V  ° 


(2.10) 


where  V  and  V  are  the  unit  cell  volumes  and  Xfl  and  are 

3  D 

the  measured  susceptibilities  of  the  III-V  analogs. 

For  the  relative  birefringence  we  obtain  from  Eq.  (2.9)  that 


An  n  -  n 

e  0 


X  -  xn 
e  0 

2(1  +  X  ) 
'  av' 


7a(T  +  a)  +  7b(T  -  o) 


6(afl  +  ab 


(2.11) 


since  X  is  always  much  greater  than  one.  For  the  II-IV-V-  compounds , 
av 

t  and  a  have  about  the  same  magnitude  so  we  can  approximate  Eq.  (2.11) 


1  7 

'a 

6a  +a 

a  d 


(t  +  a) 


(2.12) 


This  gives 


According  to  Eq.  (2.12)  the  birefringence  depends  mainly  on  the  aniso¬ 
tropic  polarizability  of  the  AC  bond.  Consequently,  since  A  is  the 
cation  with  the  largest  radius,  the  birefringence  of,  for  example, 
CdGeP2  and  CdSiP2  should  be  determined  mainly  by  the  anisotropic  polar¬ 
izability  of  the  CdP  bond.  However,  it  is  not  possible  to  estimate 
the  birefringence  of  CdSiP2  from  the  birefringence  CdGePg  because  the 
anisotropic  polarizability  of  the  CdP  bond  changes  when  silicon  re¬ 
places  germanium.  Therefore,  in  practice  the  Eqs.  (2.12)  and  (2.13) 


cannot  be  used  to  estimate  the  birefringence. 


For  reference,  Table  1  gives  a  listing  of  the  optical  properties  of 
the  chalcopyrites .  The  numbers  in  parentheses  have  not  been  measured, 
but  we  believe  they  represent  good  estimates.  The  two  most  useful 
crystals  for  nonlinear  optics  will  probably  be  CdGeASg  and  ZnGePg 
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CHAPTER  III 


BANDSTRUCTURE  OF  THE  CHALCOPYRITES 
IN  THE  VICINITY  OF  k  =  0 

CdGePg  with  a  bandgap  energy  of  I.72  eV  is  expected  to  be  the 

11-IV-V.j  chalcopyrite  compound  with  the  largest  direct  energy  gap.^ 

All  the  chalcopyrites  with  smaller  energy  gaps,  however,  are  direct 

bandgap  material  with  the  bandgap  located  at  k  =  0  .  Rowe  and 
17 

Shay  discuss  the  splitting  and  symmetry  of  the  valence  and  con¬ 
duction  bandstand  they  model  the  chalcopyrite  structure  as  a 
strained  version  of  the  III-V  sphalerite  structure.  This  gives 
good  agreement  with  the  observed  electroreflectance  spectra.  The 
S-like  conduction  band  is  singly  degenerate  (doubly  counting  spins). 
The  triply  degenerate  P-like  valence  band  (sixfold  counting  spins) 
is  split  into  three  sub-bands  because  of  spin-orbit  and  crystal 
field  interaction. 

Knowledge  of  the  splitting  and  the  curvature  of  the  valence 
and  conduction  bands  near  k  =  0  is  very  important  for  the  under¬ 
standing  of  optical  and  electrical  transport  properties.  From  the 
curvature  of  the  bands  we  determine  the  effective  masses  which 
allow  us  to  calculate  the  intrinsic  carrier  density,  the  free  »  .ectron 
absorption  cross  section,  and  the  plasma  frequency.  From  the  plasma 
frequency  we  can  determine  how  much  the  presence  of  free  carriers 
perturb  the  indices  of  refraction.  The  components  of  the  mobility 
and  conductivity  tensor  are  related  to  the  anisotropy  of  the  effective 
masses . 
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In  p-type  material  there  is  optical  absorption  due  to  introband 
transitions  between  the  valence  bands.  Our  band  structure  calcu1ation 
allows  us  to  estimate  the  magnitude  of  the  absorption  cross  section. 

In  Chapter  V  we  also  discuss  third  order  nonlinear  processes  in  the 
chalcopyrites .  The  largest  contribution  to  the  third  order  sus¬ 
ceptibility  tensor  for  III-V  compounds  with  a  free  electron  con- 
centration  of  10  cm  ^  is  from  the  nonparabolic ity  of  the  conduction 
band  near  k  =  0  similar  calculations  are  also  possible  for 

the  chalcopyrites  except  they  are  more  complicated,  since  the  lack 
of  cubic  symmetry  leads  to  much  more  complex  bandstructure.  It  is 
interesting,  however,  that  measuring  the  third  order  susceptibility 
tensor  components  provides  a  check  on  the  bandstructure  theory. 

r 

The  bandstructure  of  the  III-V  compounds  near  the  Brillouin 

zone  center  is  very  well  described  by  a  semi-empirical  theory  due  to 

2o  5 

Kane.  Borshchevskii  et  al .  adopt  these  results  directly  to  the 

chalcopyrites  to  calculate  the  effective  masses.  This  is  not  very 

satisfactory  since  there  is  a  significant  crystal  field  splitting 

21 

in  the  chalcopyrites.  We  will  therefore  extend  Kane's  theory  to 
the  chalcopyrites  and  calculate  the  energy  and  wave  functions  to 
second  order  in  k 

Including  spin-orbit  coupling,  the  Schrtidinger  equation  for  .on 
electron  in  a  periodic  potential  V(r)  is 

r  p2  *  l 

—  +  v  +  g  2  X  p),cr  ‘^k  =  Ek\  > 

.  2m  b  c 
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where  p  is  the  momentum  operator  and  <y  is  the  spin  operator. 

Introducting  the  Bloch-f unction  \k  =  u  e  where  u  has  the 

K  k  k 

periodicity  of  the  lattice,  we  obtain 


2 

1>  ft  _  _ 

—  +  V  +  -  k  •  p  f 

.  2m  m 


2  2  (^  *  p)*0-  +  p-p  (W  X  k )  * cr  |  u. 

n  c  4m  c  I  k 


(3.2) 


2  2  2  —  — 

We  neglect  the  term  (ft  /hm  c  )(W  X  k)«o-  u^  ,  since  it  is  very 
small  according  to  Kane,  because  most  of  the  spin-orbit  interaction 
occurs  in  the  interior  of  the  atom  where  the  atomic  momentum  p  is 
much  larger  than  the  crystal  momentum  fik  .  We  write  the  potential 


as  ^  =  ^0  +  Vcr  '  Here  Vq  has  cubic  symmetry  and  V 
potential  caused  by  the  tetragonal  distortion.  We  let 


if  the 

(3-3) 


be  the  Hamiltonian  for  the  unperturbed  problem  and  treat 


^  ft 

,  vcr  +  -  k-f  +  -g-g  (VV  x7).^  (5.1,) 

m  4m  c 

as  the  perturbation  where  the  tetragonal  potential  is  neglected  in 
the  spin-orbit  coupling  term.  As  basis  functions  we  take  the 
wavefunctions  for  the  valence  and  conduction  band.  Restriction  to 
these  eight  wave  functions  means  that  the  influence  of  the  higher 


and  lower  lying  band*  are  considered  call.  This  approximation  works 

bast  when  the  energy  separation  to  the  other  bands  at  k  =  0  is 

"uch  larger  than  the  bandgap  energy.  We  follow  Kane  and  use  the 
basis  functions: 

1 iSi)  ’  l£l/2  (*  -  »)f>  ,  |21>  ,  I  2'1/2(X  +  iY)r>  , 

I ‘St)  ,  I21/2  (X  +  1Y)1)  ,  Izt)  ,  |  2-‘/2(x  .  iy)i) 

four  functions  are  respectively  degenerate  with  the  last 
four.  The  symbols  T  and  i  indicate  spin  up  and  down  referred  to  the 
n-axls.  s  refers  to  the  conduction  band  wave  functions.  They  trans¬ 
form  as  atomic  s-functlons  under  the  symmetry  operations  of  the 
tetrahedral  group.  Similarly  X  ,  y  ,  and  z  refer  to  the  valence 
band  wave  functions,  and  they  have  the  symmetry  properties  of  the 

atomic  p-functions  *  ,  y  ,  and  r  under  the  tetrahedral  group 
operations. 

with  these  basis  functions,  it  follows  that  the  8  X  8  Hamil- 
tonian  matrix  can  be  written 


(5-5) 
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where 


E  ♦  *-$ 
P  3  3 


~A 


21/2(k  +  ik  )P 

x  y 


2l/2(kx  -  iky)P 


2 

E  -  f  6 
P  3 


E  +  I  +  f 

P  3  3 


(3.6) 


■2  /2(kx  +  iky)p  0  2  7  (kx  -  iky)p 


(3-7) 


Here  E  and  E  are  the  eigenvalues  of  the  unperturbed  Hamiltonian, 
s  p 

referring  to  respectively  the  two-fold  degenerate  conduction  band 
and  the  six- fold  degenerate  valence  band.  We  further  have  defined 
the  quantities 


P  =  -  i  -  <S|pJZ> 
m 
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(3.8) 


with 


(S|pz|Z)  =  <s|pjx>  =  <s|py|Y>  ; 


5»i  av 

a  =  (x|  —  P  -  — 1  p  |Y> 

4m  c“"  Sx  y  ay 


(3.9) 


with 


av.  av.  avA  av 

-  (x|r^  -^|2> 

av  avA 

=  |  — 2-  p  — —  p  |z> 

ay  z  az  y 


and 


<x|v 

'  1  cr 

|x>  = 

<S|Vct 

|S> 

♦i 

<Z  I V 
'  1  cr 

jz)  = 

<S|Vcr 

|S> 

2& 

3 

(3.10) 


We  have  not  included  the  term  <S|V  jS)  in  Eq.  (3.6)  since  it  only 
shifts  the  eigenvalues  by  a  constant.  To  determine  the  level  splitting 
near  k  *  0  to  first  order  in  H'  because  of  the  spin-orbit  coupling 
A  and  the  crystal  field  splitting  &  ,  wa  diagonalize  the  matrix 

in  Eq.  (3-5) •  We  find  that  the  conduction  band  remains  two-fold 
degenerate  while  the  valence  band  splits  into  three  doubly  degenerate 
levels.  The  four  doubly  degenerate  eigenvalues  are  determined  by  the 
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secular  equation: 


2A&  ' 

2  / 

E 

('  -  V 

E2  +  (6  +  A)E  +  - 

3 

-  <v>  (E  +  -j 

2 


0  , 


(3.11) 


where  E  =  Ek  -  (ft2k2/2m)  .  We  have  fixed  the  zero  on  the  energy 
scale  in  Eq.  (3-H),  bY  letting  Ep  +  &/3  +  A/3  =  0  .  At  k  =  0 
the  eigenvalues  are  Eg  >  Ei  >  0  >  and  E2  Wit^ 

E1  2  =  -  |  (A  +  &)  1  \  \j(b  +  A^2  ‘  ^  *  ^‘l2) 

In  the  chalcopyrites  we  have  A  >  0  and  6  <  O17  which  that 

E^  >  0  and  Eg  <  0  .  E^  is  therefore  the  energy  of  the  highest 

lying  valence  band  and  the  bandgap  energy  is  EG  =  Es  "  El  *  We 
use  this  together  with  Eq.  (3-12)  to  rewrite  Eq.  (3-H).  We  find 

that 

E(E  -  Ej_  -  Eg)(E  -  Ex) (E  -  Eg)  -  i**)2*  ^E  + 

+  (k2  +  k2)P2  [(E  -  EX)(E  -  Eg)  -  |(E  +  & 

(3-13) 

For  small  k2  ,  the  solutions  of  Eq.  (3*13)  give  parabolic  energy 
bands.  Including  the  (h2k2/2m)  term  from  Eq.  (3-2),  we  obtain 
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to  second  order  in  k  that  the  energies  of  the  conduction  and  valence 
bands  are 


E 
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(MM 

The  above  equations  give  the  energy  splittings  and  the  effective  masses 

of  the  four  bands  excluding  the  effect  of  higher  lying  bands.  It  is 

worth  noting  that  the  expressions  for  Ey  and  Ey  are  not  valid 

in  the  limit  6  -*  0  .  We  now  proceed  to  calculate  the  wave  functions. 

The  lattice  periodic  functions  |u  )  are  given  by  a  linear  combination 

m 

of  the  eight  basis  functions  |n)  .We  have 

IV  ’  I  aJn>  <  (5.15) 

n 
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and  we  determine  the  coefficients  a^  by  solving  the  eight  homo¬ 
geneous  equations 

V  a  (H  -  E  6  )  =  0  ;  s  =  1  to  8  (3.16) 

mnl  sn  m  sn j  *  ' 

n  ' 

for  each  eigenvalue  E  of  the  Hamiltonian  matrix  H  in  Eq.  (3.5I. 

m  sn 

We  then  obtain  Eqs.  (3-17)  and  (3.I8) 
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In  Eqs .  (3.I7)  and  (3-18)  A  and  B  distinguish  between  the 
two  wave  functions  of  each  doubly  degenerate  band.  The  Eqs.  (3.14), 
(3.17),  and  (3. 18)  are  the  main  results  of  this  chapter.  With  these 
equations  we  can  calculate  energy  splittings,  effective  masses,  and 
optical  transition  probabilities  between  the  bonds.  The  equations 
depend  only  on  the  following  four  parameters:  the  bandgap  energy 
Eg  ,  the  spin-orbit  splitting  energy  A  ,  the  crystal  field  splitting 
parameter  5  ,  and  finally  the  matrix  element  P  .  For  the  band- 

gap  energy  we  use  the  measured  value.  The  spin-orbit  coupling  and 
crystal  field  splitting  energies  have  been  measured  for  four  of  the 

II- IV-V2  compounds:  CdGeP2,  CdSnP2,  ZnSiAs2,  and  CdSiAs^16  For 
the  other  chalcopyrites  we  estimate  the  spin-orbit  coupling  from  tne 

III- V  compounds  by  averaging  the  experimental  values  for  the  III-V 

22 

analogs.  For  the  crystal  field  splitting  we  use  the  results  of 
a  pseudopotential  calculation, 21  noting  that  with  our  definition  of 
the  crystal  field  splitting  parameter  it  assumes  negative  values 
in  the  chalcopyrites.  It  is  also  possible  to  estimate  the  crystal 
field  splitting  by  considering  the  chalcopyrites  as  stressed  versions 
of  the  III-V  compounds  as  suggested  by  Shay  et  al.1^'2^  A  comparison 

pli 

with  Poliak  et  al.  yields 

3  ,  5 

6  =  ~~bE001  =  -  5b^zz  ■  exx>  =  ”bT  >  (3.19) 

where  b  is  the  deformation  potential  of  the  III-V  analogs  for  stress 
in  the  [001]  direction.  For  GaAs  b  =  -  1.75  eV.25,  Since  we  do  not 
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know  b  for  InAs ,  we  use  this  value  to  estimate  &  =  -  0.29  eV  for 
CdGeAs^  which  is  close  to  -0.25  predicted  by  the  pseudopotential 
calculation. 

The  matrix  element  P  is  nearly  constant  for  all  the  III-V 
compounds  that  can  be  described  by  Kane's  theory.  Within  20$  we 
have 


2  * 

P  =  —  X  20  eV 
2m 


(3.20) 


for  the  III-V  compounds,  and  we  assume  this  value  for  the  chalcopyrites , 
Close  to  k  =  0  the  enbrgy  bands  are  of  the  form 


ft 

E  =  E0±- 


k?  +  k2  k2 
x  y  z 

- u-  +  - 


(3.21) 


"T  "L  J 

with  the  plus  and  minus  sign  referring  to  respectively  the  conduction 

and  valence  bands.  In  Eq.  (3.21)  itVj,  and  m^  define  the  transverse 

and  tfcsu  longitudinal  effective  masses.  We  can  use  these  to  define 

a  denfity  of  states  effective  mass  m, 

de 


nr 


de 


2 

mLmT 


(3-22) 


(3-23) 


and  a  conductivity  effective  mass  m 

a 

1  1  / 1  2 

"V  5'\”l  "l, 

(cf.  reference  27).  Table  II  gives  the  calculated  effective  masses 

for  CdGeAs  and  CdGeP  together  with  the  values  for  E„  A  and  & 

<-  G  9 
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TABLE  II 


CALCULATED  VALENCE  BAND  SPLITTING  ENERGIES  AND  EFFECTIVE 
MASSES  FOR  THE  CONDUCTION  AND  VALENCE  BANDS 


CdGeASg 

CdGeP2 

eg 

— 

0.53 

1.72 

A 

0.38 

0.11 

8 

-  0.25 

-  0.20 

El 

0.20 

0.17 

E2 

-  0.32 

-  0.08 

mT 

mL 

m, 

de 

“tT 

“l 

m, 

de 

m<r 

c 

0.039 

0.028 

0.035 

0.034 

0.088 

0.079 

0.085 

0.085 

V1 

0.77 

0.031 

0.26 

0.087 

- 

0.097 

v2 

0.079 

— 

0.23 

— 

V3 

0.14 

0.69 

0.24 

0.19 

0.24 

All  energies  are  in  eV  and  the  effective  masses  are  in  units  of  the 


free  electron  mass. 


used  in  the  calculation.  The  table  also  lists  the  splitting  energits 
Ei  and  of  the  valence  band.  We  have  assumed  direct  bandgaps, 

setting  Eg  equal  to  the  measured  bandgap  energy.  There  is,  however, 
some,  uncertainty  concerning  CdGeAsg.  Early  reports  of  a  direct  band- 
gap'  P  have  been  questioned  recently,29  although  the  fact  that  both 
of  the  III-V  analogs  GaAs  and  InAs  have  direct;  oandgaps  tends  to 
support  the  notion  of  a  direct  bandgap.  With  this  qualification  in 
mind.  Fig.  3  shows  the  expected  bandstructure  of  CdGeAs^  near  k  =  0. 
rhe  highest  valence  band  has  a  smaller  longitudinal,  but  a  larger 
transverse  effective  mass  than  the  v2-band.  A  few  of  the  effective 
masses  of  the  valence  bands  are  not  listed  in  Table  II.  These  can 
only  be  determined  accurately  by  considering  the  effect  of  higher 
bands  since  they  receive  very  little  contribution  from  the  conduction 
band  interaction.  No  measured  effective  masses  have  been  reported 
ror  CdGePg.  In  CdGeAs2,  however,  infrared  plasma  reflection50 
and  thermoelectric  power  measurements28'51  have  been  used  to  determine 
the  electron  mass.  The  effective  mass  depends  on  the  carrier 
concentration.  Extrapolation  to  small  carrier  concentrations  results 
in  an  effective  mass  of  0.02  to  0.03  times  the  free  electron  mass. 

This  is  in  reasonably  good  agreement  with  the  calculated  effective 
mass  of  0.035  in  Table  II.  There  is  only  one  reported  measurement 
of  the  effective  hole  mass  in  CdGeAs,,.51  The  mass  was  deduced  from 
measurements  of  the  thermoelectric  power.  This  gave  a  density  of 
state  effective  mass  of  0.3  which  agrees  surprisingly  well  with  our 


calculated  value  of  0.26.  We  therefore  conclude  that  the  generalized 
Kane  model  that  we  have  described  gives  adequate  description  of 
the  bandstructure  near  k  -  0  in  CdGeAs2  and  probably  works  well  also 
for  the  other  chalcopyrites . 

Finally  in  this  chapter,  we  list  the  matrix  elements  for  optical 
transitions  between  the  highest  and  the  two  lower  valence  bands. 

In  Chapter  VI  we  will  refer  back  to  the  results  when  we  estimate 
the  magnitude  of  the  absorption  in  p-tvpe  material.  (See  Eq.  (3.24) 
on  the  following  page). 
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CHAPTER  IV 


ELECTRIC  PROPERTIES  AND  LINEAR  EXPANSION  COEFFICIENT  OF  CdGeAs, 

Several  absorption  mechanisms  in  the  chalcopyrites  are  related 
to  the  semiconductor  properties.  Therefore  knowledge  of  the  electric 
transport  properties  is  necessary  for  a  full  evaluation  of  the  material. 
It  is  well  known  that  controlled  doping  may  improve  the  optical 
quality  of  semiconductors.  As  an  example,  doping  GaAs  by  chromium 
introduces  deep  acceptor  levels  that  remove  all  the  free  electrons. 

This  results  in  high  resistivity  material  of  good  optical  quality. 
Another  example  is  CdSe  where  Se  compensation  reduces  the  carrier 
concentration.  Besides  causing  scattering  leading  to  absorption,  the 
presence  of  free  carriers  perturbs  the  indices  of  refraction 
and  thereby  the  phasematching  angle  for  nonlinear  optical  inter¬ 
actions.  An  additional  absorption  mechanism  exists  in  p-type 
material  due  to  intraband  transitions  between  the  split  valence 
bands.  This  absorption  mechanism  has  a  magnitude  proportional  to  the 
free  hole  concentration.  Another  common  cause  of  optical  absorption 
in  semiconductors  is  impurity  absorption.  From  the  temperature 
dependence  of  the  Hall  coefficient  we  can  derive  information  about 
the  position  and  concentration  of  the  impurity  levels,  and  by  combining 
the  Hall  coefficient  with  the  measured  conductivity  we  can  calculate 
the  carrier  mobility.  The  mobility  is  sensitive  to  the  general 
crystal  quality  and  increases  with  a  smaller  number  of  crystal 
imperfections. 
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In  this  chapter  we  discuss  the  electrical  transport  equations 
necessary  to  interpret  the  experimental  results  and  then  present 
the  resuits  for  p-type  CdGeAs^.  We  will  also  discuss  the  measure¬ 
ments  of  the  linear  thermal  expansion  coefficient  for  CdGeAs^.  The 
expansion  coefficient  is  strongly  anisotropic  and  explains  at  least 
partially  the  crystal  cracking  during  growth. 

A.  ELECTRIC  TRANSPORT  PROPERTIES 

1 .  Theory 

Introducing  an  electron  mobility  ^  and  a  hole  mobility  p  , 
we  have  that  the  electric  conductivity  is  given  by 

cr  =  e(n^c  +  p^)  ,  (l*.l) 

when  there  are  n  electrons  in  the  conduction  band  and  p  holes  in 
the  valence  band.  The  mobility  depends  on  the  collision  time  of  the 
carriers  and  the  effective  mass*  For  holes  we  have 

e  (t) 

=  -  >  (fc.2) 

m 

v 

and  a  similar  expression  holds  for  the  electrons,  (t)  is  the 
collison  time  averaged  over  the  carrier  velocity  since  t  in 
general  is  velocity  dependent. 

In  the  case  of  an  anisotropic  valence  or  conduction  band,  <7  is 
a  tensor.  From  band  structure  calculations  near  the  Brillouin 
zone  center  we  have  shown  that  the  chalcopyrites  have  a  nearly 
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isotropic  conduction  band,  but  that  the  valence  band  is  strongly 
anisotropic.  It  is  convenient  to  define  a  longitudinal  and  a 
transverse  conductivity  for  conduction  parallel  and  perpendicular  to 
the  c-axis  direction.  For  p-type  material  we  have  according  to 
Eqs.  (l+.l)  and  (1+.2)  that 


(M) 


where  refers  to  the  top  valence  band.  Using  the  effective  masses 

in  Table  II,  we  obtain  a =  25  for  CdGeAs^  The  ratio  between 
the  electron  and  hole  mobility  depends  on  which  scattering  mechanism 
limits  the  collision  time.  The  average  collision  time  is  in  general 
different  for  electrons  and  holes.  Pure  acoustic  scattering  is  most 
important  at  high  temperatures  due  to  the  T_3/2  temperature  dependence. 

In  that  case  the  ratio  between  the  longitudinal  components  of  the  mobility 

?7 

tensor  is  1 


and  we  obtain  a  similar  expression  for  the  transverse  components. 

The  subscript  de  in  Eq.  (k.k)  refers  to  the  density  of  states  effective 

mass.  Applying  Eq.  (1».4)  to  CdGeAs,,  we  obtain  ^  =  22  and 

Mc  T^v  =  •  *  1,L 

>  V1,T 
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By  a  slight  extension  of  the  results  in  reference  27,  we  obtain 
expressions  for  the  Hall  constant  R  in  crystals  such  as  CdGeAs^ 
with  spheroidal  constant  energy  surfaces  and  a  direct  bandgap  at 
k  =  0  .  For  the  magnetic  field  applied  along  the  c-axis  we  find 


2  2 

-L  PMvJ  "  n>1c,T 

|e|  (puV)I  +  n nc>T)‘ 


Ct-5) 


and  for  the  magnetic  field  perpendicular  to  the  c-axis  we  have 


V  =  , 


r  PMv  T^v  L  "  n^c  T^c  L 

|e|  (PHVjT  +  np^T)(pp^L  +  np^L) 


(4.6) 


2  2 

In  the  above  equations  r  =  (t  )/(t)  is  a  constant  close  to  unity, 

i.e.  r  =  1.18  for  pure  acoustic  scattering  and  r  =  1-93  for 

27 

ionized  impurity  scattering.  With  only  one  type  of  carrier  the 
equations  simplify.  For  p-type  material  we  have 


(4.7) 


and  in  this  case  the  Hall  constant  is  a  dir  .ct  measure  of  the  carrier 
concentration. 

In  intrinsic  semiconductors  the  number  of  holes  equals  the  number 
of  electrons.  With  the  definitions 


Nc  =  2 


2nkTm 


c  ,de  I  2 


(4.8) 
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N 


v 


2 


27TkTm 


c  ,de 


1 

2 


(h.9) 


we  calculate  the  intrinsic  carrier  concentration  nt  by  the  equation 


n'T  =  N  N  exp  ( - - 

1  cv  '  kT, 


(4.10) 


Using  the  effective  masses  mc  =  0.035  m  and  mv  jg  “  0.26  m  , 
we  obtain  for  CdGeAs0  that 

Figure  1*  shows  the  temperature  dependence  of  the  intrinsic  carrier  con 

centration.  At  300°K  we  have  N  =  1.64  x  10^  cm  ^  and 

c 

18  -7j 

N  =  3-32  X  10  cm  J  . 
vi 

For  extrinsic  semiconductors  assuming  nondegenerate  bands,  the 
product  np  of  the  electron  and  hole  concentration  is  constant 
independent  of  the  impurity  concentration  and  it  is  given  by 


=  1.42  X 


I0lk  T 


3/2 


exp  - 


np  = 


(4.12) 


The  semiconductor  is  nondegenerate  when  the  Fermi  level  is  in  the  for¬ 
bidden  energy  gap  and  is  separated  from  the  valence  and  the  conduction 
band  by  an  energy  of  more  than  kT  .  For  p-type  material  this 
condition  implies  that  p/Nv  <  l/e  with  e  =  2.718  .  This  condition 
can  be  rewritten  such  that  for  a  carrier  concentration  p  the  semiconductor 


is  nondegenerate  at  temperatures  larger  than  the  degeneracy  temperature 


Tp  defined  by  the  equation 


td  =  ^ 


2nk  m 


(4.13) 
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For  a  hole  concentration  of  10  cm  J  in  CdGeAg2  the  degeneracy  tem¬ 


perature  is  as  low  as  12°K 


Finally  in  this  section,  we  discuss  the  temperature  variation 


of  the  carrier  concentration  in  a  material  with  N,  donors  and  two 

d 


acceptor  levels  of  concentration  N  and  N  respectively.  The 


results  apply  to  the  measured  p-type  CdGeAs2  samples.  Tte  charge 


balance  leads  to  the  equation 


n  +  '  ",  )  -  P  +  (N,  -  ".)  , 


(4.U) 


where  nd  is  the  number  of  filled  donor  levels  (un-ionized  donors) 


and  n„  and  are  the  number  of  unfilled  acceptor  levels  (un¬ 


ionized  acceptors).  We  consider  p-typa  naterial.  The  donors  are 


ionized  such  that  nd  =  0  .  Furthermore,  we  let  refer  to 


a  shallow  acceptor  level  also  fully  ionized  (nn  =  0)  .  N  refers 

1  a2 


to  a  deeper  partially  ionized  acceptor  level.  For  an  acceptor  ioniza¬ 


tion  energy  of  £  ,  the  number  of  un-ionized  acceptors  is  given  by 


=  - - C _ 

1  +  |  exp  [(Ej,  -  e  )/kT] 


(*.15) 


■ 

[  ' 


' 


Here  the  Fermi  energy  is  measured  from  the  top  of  the  valence  band  and 
the  factor  of  two  is  the  spin  degeneracy.  For  simplicity  we  also 
neglect  n  in  Eq.  (4.14).  This  is  allowed  as  long  as  the  hole 
concentration  is  much  larger  than  the  intrinsic  carrier  density. 

For  a  nondegenerate  semiconductor,  the  concentration  of  holes  in  the 
valence  band  is  given  by 

-FT/kT 

P  =  V  ,  (4.1.6) 

which  substituted  into  Eq.  (4.15)  yields 


2  1  +  2  (Nv/p)  exp  (“  ea  /kT) 


(4.17) 


This  eliminates  the  Fermi  energy  from  Eq.  (4.14)  and  we  can  solve  for 

the  hole  concentration.  We  obtain 
i  .  -e_  /kT 

p  -  a  <%  -  Kd  -  5  V  8  ) 


+  2  -  "d  -  I  V  W  +  2(Mai  +  ^  -  »d)«v.'eaa/kT  I5 

(4.16) 

with  the  asymptotic  solutions  p  =  N  -  Nd  at  low  temperatures  and 
P  =  Nfli  +  N&2  -  Nd  at  high  temperatures.  The  donors  partially 
compensate  the  hole  concentration.  At  intermediate  teinperatures 
there  is  an  approximate  solution  provided  the  condition 


ea >T  I  5' 


Na  "  Nd  «  P  «  Na  +  N  -  N . 
al  d  ai  a2  d 


(U. 19) 
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can  be  satisfied  which  is  only  possible  when  N  »  N  -  N 

a2  ai  d  ' 

With  Eq.  (4.19)  satisfied,  the  approximate  carrier  concentration  is 
given  by 


%  N  N 
2  v 


-£  /2kT 
s  a2 


(4.20) 


2.  Experiment 


We  have  measured  the  resistivity  and  Hall  constant  between  77°K 
and  400°K  for  p-type  CdGeAs,,.  The  samples  came  from  boules  grown  from 
stoichiometric  melt.  For  most  samples  we  used  the  van  der  Pauw 
m.thod,52,33  We  did  not  orient  the  samples.  They  consisted  of  one  or 
a  few  crystallites  and  the  sample  diameter  to  thickness  ratio  was 
about  10.  It  would  have  been  desirable  to  use  oriented  rectangular 
bar  shaped  samples  since  this  would  have  allowed  us  to  measure 
the  individual  components  of  the  conductivity  tensor.  However, 
due  to  the  extensive  cracking  it  was  not  possible  to  obtain  suffi¬ 
ciently  large  oriented  single  crystal  samples.  We  used  indium 
to  form  ohmic  contacts.  The  contacts  worked  well  down  to  liquid 
nitrogen.  Small  balls  of  indium  were  squeezed  onto  the  sample  by 
a  teflon  coated  tv-er.er  and  then  alloyed  in  an  Hg  atmosphere  for  5 
to  10  minutes  at  350 °C. 

Figure  5  shows  the  temperature  dependence  of  the  Hall  constant. 

We  used  a  current  of  0.5  mA  and  a  magnetic  field  of  2000  Gauss. 

This  gave  a  Hall  voltage  of  about  100  ^V.  The  Hall  constant  is 
almost  independent  of  temperature  between  77°K  and  250°K.  At 
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FIG.  5— Temperature  dependence  of  the  Hall  constant  for  p-type  CdGeAs 


higher  temperatures  it  starts  to  decrease.  The  resistivity  plotted 
in  Fig.  6  also  has  the  same  general  temperature  dependence.  It  is 
constant  between  77°K  and  250°K  and  decreases  at  higher  temperatures. 
This  agrees  with  published  data.^  ***  Vaipolin  et  al.  observed 
that  the  Hall  constant  started  to  decrease  above  200°K  and  claimed 
that  this  was  due  to  the  onset  of  mixed  conductivity.  This,  however, 
is  not  compatible  with  their  results  showing  the  Hall  constant  goes 
to  zero  at  500°K.  With  a  bandgap  of  0.55  eV,  it  is  easy  to  estimate 
from  Eqs .  (4.5)  or  (4.6)  and  Fig.  4  that  the  temperature  difference, 
measured  in  units  of  1000/T,  between  the  onset  of  mixed  conductivity 
and  zero  Hall  constant  should  be  less  than  0.5°K  1.  This  is  much 
smaller  than  the  difference  between  200°K  and  500°K  which  is  3  K  . 
Furthermore,  the  resistivity  changes  too  quickly  with  the  temperature 
to  be  explained  by  purely  a  change  in  the  mobility.  The  correct 
interpretation  is  probably  that  the  crystals  have  two  acceptor  levels, 
one  shallow  and  one  deep.  The  shallow  level  is  completely  ionized 
above  77°K.  Between  77°K  and  250°K  we  have  a  region  of  saturated 
conduction  with  the  carrier  concentration  equal  to  -  Nd  . 

Above  250°K  the  effect  of  the  deep  acceptor  level  starts  to  be  notice¬ 
able. 

When  there  is  only  one  type  of  carrier,  the  Hall  constant  gives 
a  direct  determination  of  the  carrier  concentration.  Neglecting  the 

constant  r  ,  it  follows  from  Eq.  (4.7)  that  p  =  l/|e|R  .  This 

o  15  **3 

gives  a  carrier  concentration  at  77  K  of  6.6,  8.1,  and  5*3  *  40  cm 
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FIG.  6 — Temperature  dependence  of  the  resistivity  for  p-type  CdGeAs 


respectively  for  the  samples  from  boules  b6,  BJf  and  b8.  Since  the 

samples  have  a  nondegenerate  carrier  concentration,  we  can  use  £q. 

(4.20)  to  obtain  a  rough  estimate  of  the  impurity  concentration  that 

gives  rise  to  the  deep  acceptor  level.  Equation  (4.20)  gives  the 

carrier  concentration  at  intermediate  temperatures  and  we  determine 

the  ionization  energy  of  the  deep  levels  from  the  slope  in  Fig.  5* 

This  gives  e  «  0.2  eV,  and  we  est 'mate  the  impurity  concentration 

a2 

17  -*5 

to  be  Na,.j  «  4  X  10  cm  J  for  sample  B8. 

The  Hall  constant  and  the  resistivity  allow  us  to  determine  the 
Hall  mobility  which  is  defined  by  ^  =  R/p  •  For  our  samples  the 
Hall  mobility  is  temperature  independent  up  to  almost  room  tempera¬ 
ture.  It  then  starts  decreasing.  The  Hall  mobilities  calculated  at 
77°K  for  samples  B6,  B7,  and  B8  are  respectively  676,  606,  and  1540 

p 

cm  /Vsec.  The  difference  between  B8  and  the  two  other  samples  might 

be  partially  due  to  orientation  effects.  For  comparison  we  also  list 

^8 

the  previously  published  results.  Reference"  gives  a  hole  mobility 

p 

of  20  cm /Vsec  and  reference  34  lists  the  electron  and  hole  mobilities 

p  o 

as  1000  cm  /Vsec  and  240  cm  /Vsec  respectively.  We  therefore  conclude 
that  our  material  is  probably  of  slightly  higher  electronic  quality 
with  fewer  ionized  impurities  acting  as  scattering  centers  to  limit 
the  mobility. 

B .  MEASUREMENT  OF  THE  LINEAR  THERMAL  EXPANSION  COEFFICIENT 

We  have  measured  the  linear  thermal  expansion  coefficient  for 
CdGeAs g  in  the  temperature  range  80°C  to  170°C,  using  an  optical 
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interference  technique.  Figure  7  shows  the  experimental  setup.  To 
avoid  vibrations  the  experiment  should  be  performed  on  a  stable  table. 

A  small  tilt  between  the  microscope  slides  gives  rise  to  fringes, 
and  the  fringe  spacing  on  the  screen  decreases  with  increasing  tilt 
angle.  The  microscope  slides  were  coated  to  obtain  good  fringe 
contrast.  At  times  the  contrast  could  be  improved  by  a  wedge 
blocking  the  zero  order  reflection.  We  had  regular  air  in  the  furnace. 
Above  200°C  it  is  necessary  to  use  an  inert  atmosphere  or  vacuum  to 
avoid  growth  of  an  oxide  layer.  A  vacuum  has  the  additonal  advantage 
that  it  is  otherwise  necessary  to  correct  for  the  change  in  the  index 
of  refraction  with  the  temperature  of  the  gas  between  the  microscope 
slides . 

Let  t  be  the  sample  thickness  and  n  the  index  of  refraction 
for  air.  The  phase  difference  associated  with  one  double  reflection 
between  the  microscope  slides  is  then  given  by  cp  =  (lm/x)ni 
Here  X  =  6328  X  is  the  wavelength  of  the  He-Ne  laser.  A  temperature 


We  measure  Z^cn  by  observing  the  fringe  movement  on  the  screen.  For 
s  fringes  passing  through  a  reference  point,  we  have  =  27TS  .  By 
substituting  this  into  Eq.  (4.21)  we  determine  the  linear  thermal 


coefficient  a  .  We  find 


20  cm  f  *  1.45  cm 


FIG.  T — Experimental  setup  for  measuring  the  linear  thermal  expansion  coefficient. 


We  see  from  Eq .  (4.22)  that  we  need  to  know  dn/dT  for  air.  The  index 
of  refraction  for  air  can  be  written  as  n  =  1  +  ^  X  with  the  sus¬ 
ceptibility  X  proportional  to  the  density  N/v  =  p/kT  .  From 
this  it  follows  that 


(4.23) 


For  p  =  Pq  and  defining  3  =  1/Tq  ,  we  can  rewrite  this  equation  as 

n_  -  1 
T 

— -  =  1  +  p(T  -  TQ)  .  (4.24) 

°T  '  1 


Letting  =  273°K  ,  we  should  theoretically  have  3  =  0.00367 

Experimentally  3  is  found  to  be  equal  to  0.003679  fcr  X  =  6328  8 

o  35 

and  Tq  =  273  K.  For  the  same  wavelength  and  temperature  and  a 

oressure  of  760  mm  Hg  the  index  01  refraction  for  dry  air  is  n_  = 

r0 

1.0002921  .  The  temperature  dependence  is  obtained  by  differentiating 

Eq.  (4.24).  We  obtain 


T  T 

-  =  -  p  -  .  (4.25) 

dT  1  +  3(T  -  T  ) 


The  calculated  dr^/dT  is  plotted  in  Fig.  8. 

In  order  to  test  the  method  we  first  measured  a  crystal  with  a 
known  expansion  coefficient.  We  used  a  single  crystal  of  silicon 
3.99  nan  thick.  The  results  are  listed  in  Table  III.  In  Fig.  9  we 
compr *  i  the  results  with  the  data  for  polycrystalline  silicon  from 
reference  J>6.  The  agreement  is  good. 
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We  had  two  crystals  of  CdGeAs^,  one  c-axis  cut  4.01  mm  thick  and 
one  a-axis  cut  2.96  mm  thick.  CdGeASp  proved  to  have  a  very  aniso¬ 
tropic  linear  thermal  expansion  coefficient  with  Ctj^  more  than  ten 
times  larger  than  ctj|  .  The  results  for  a-axis  expansion  are  listed 
in  Table  IV.  For  the  c-axis  expansion  we  did  not  obtain  equally 
accurate  measurements  because  cr j|  was  of  the  same  order  as  dn/dT 
for  air,  but  of  opposite  sign,  such  that  the  fringe  movement  was 
almost  zero.  The  results  are  in  Table  V.  In  Fig.  10  we  plot  the 
temperature  dependence  of  the  expansion  coefficients.  The  expansion 
coefficients  increase  slightly  with  the  temperature. 


TABLE  IV 


LINEAR  THERMAL  a -AXIS 


EXPANSION  COEFFICIENT  FOR  CdGeAs 


svaOPO 


CHAPTER  V 


CRYSTAL  GROWTH  AND  MATERIAL  EVALUATION 


In  this  chapter  we  briefly  describe  the  growth  technique  and  then 
the  testing  procedures  used  to  evaluate  chalcopyrite  crystals.  The 
problems  of  main  concern  are  the  optical  transmission  and  crystal 
cracking.  Fortunately,  the  necessary  crystal  size  for  most  nonlinear 
optic  applications  is  only  a  few  millimeters  so  samples  can  be  cut  from 
the  best  sections  of  the  boules.  The  crystals  are  grown  at  the  Center 
for  Materials  Research  (CMR)  at  Stanford. 

A.  GROWTH  METHOD 
1 .  CdGeAs„ 

References  37  and  38  discuss  the  equilibrium  phase  diagram  of 
CdGeASg.  The  crystal  melts  congruently  at  670°C.  The  homogeneity 
region  is  believed  to  be  very  narrow.  It  probably  extends  only  a 
traction  of  a  percent  away  from  the  stoichiometric  composition  which 
is  helpful  in  obtaining  uniform  single  crystals.  All  our  measurements 
have  been  on  samples  grown  by  the  vertical  Bridgman-Stockbarger  tech¬ 
nique  . 

Stoichiometric  proportions  of  the  elements  (6n)  are  placed  in  a 
quartz  tube  which  then  is  sealed  under  vacuum.  The  material  is  re¬ 
acted  for  16  hours  at  a  temperature  of  730  to  7^0°C.  The  temperature 
is  raised  slowly  since  the  arsenic  vapor  pressure  is  about  10  atm  at 
730°C.  After  complete  reaction,  the  material  is  ready  for  Bridgman 
growth.  The  furnace  is  vertical  with  a  temperature  profile  as  illus¬ 
trated  in  Fig.  11.  The  quartz  crucible  is  placed  near  the  top  of  the 
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'rowth  of  CdGeAs 


furnace  and  is  lowered  slowly  through  the  temperature  gradient.  To 
aid  formation  of  a  single  crystal,  the  crucible  has  a  1  to  l|  inch  long 
capillary  at  the  bottom.  The  synthesis  and  growth  can  be  done  in  one 
or  two  steps.  The  two  step  procedure  has  the  disadvantage  that  some 
vapor  in  the  crucible  deposits  on  the  walls  when  the  synthesized 
material  is  quenched,  and  it  is  not  possible  to  retrieve  this  deposit. 
The  one  step  process  assures  that  no  material  is  lost.  However,  there 
is  a  problem  of  obtaining  uniform  mixing  and  achieving  fully  reacted 
material  in  the  crucible  as  well  as  in  the  capillary  before  the  growth 
starts.  An  ultrasonic  vibrator  aids  the  mixing  process.  To  prevent 
reaction  of  the  melt  with  the  crucible  walls  a  carbon  layer  is  used 
to  coat  the  inside  of  the  quartz  crucible.  Boules  grown  without  the 
earbon  layer  did  not  show  any  noticeable  change  in  the  optical  trans¬ 
mission. 

Figure  12  shows  a  photograph  of  a  boule.  The  boules  are  usually 
2.5  cm  long  and  1.3  cm  in  diameter  with  the  general  properties  fairly 
repeatable  from  run  to  run.  There  is,  however,  a  secondary  nuclea-* 
tion  problem.  The  boules  start  as  single  crystals,  but  after  0.5  to 
1  cm  of  growth  the  melt  becomes  unstable  and  numerous  small  crystals 
form  from  nucleation  sites.  In  trying  to  eliminate  this  problem,  the 
growth  conditions  have  been  varied  over  a  wide  range.  The  lowering 
rate  of  the  crucible  has  been  changed  between  0.25  and  27  mm/hour  and 
the  temperature  gradient  between  11.8  and  3^°C/cm,  this  has  not  resulted 
in  any  improvement.  Another  problem  is  extensive  crystal  cracking  which 
occurs  at  all  investigated  growth  conditions. 
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2-  CdGeP. 

- 2 

The  melting  point  of  CdGePg  is  at  ?90°C.  Because  of  the  high 
phosphorus  vapor  pressure,  CdP2  is  synthesized  first  and  then  reacted 
with  Ge  to  form  CdGe »g  at  a  temperature  above  the  melting  point  of  Ge. 
The  fully  reacted  mat«v ial  is  then  used  to  grow  single  crystals  by  the 
Bridgman-Stockbarger  technique  similar  to  the  growth  of  CdGeAs,,.  Due 
to  observed  Si  vitrification,  a  carbon  crucible  is  used  to  synthesize 

the  CdGeP2.  The  resulting  crystals  show  the  same  cracking  pattern  as 

observed  for  CdGeAs  . 

2 

B.  INITIAL  TESTING 
1 •  Polishing  and  Etching 

Both  CdGeAs ^  and  CdGeP2  polish  veil.  Since  the  crystals  ore 
anisotropic,  the  grain  boundaries  and  the  twin  lines  can  be  seen  on 
a  polished  surface  using  a  microscope  and  partially  crossed  polarizer.. 
Without  polarizer.  It  Is  necessary  to  use  a  suitable  etch  to  reveal  the 
structure.  A  good  etch  for  CdGeAs.,  Is  lH^:  and  etching 

for  30  sec.  p„r  CdGeP,,  etching  for  15  sec  In  a  10*  Hr,,  solution  in 
ethyl  alchol  gives  good  results.  The  best  section  of  the  boule  1, 
usually  the  region  vlthln  0.5  to  1  cm  from  the  capillary.  This  section 
most  often  consists  of  two  or  three  crystallites.  For  growth  from 

stoichiometric  melt  there  is  no  evidence  of  inclusions  except  near  the 
top  of  the  boules. 

2.  Cracking 

Due  to  the  extensive  cracking,  the  largest  single  crystals  of 
CdGeAs 2  useful  for  nonlinear  optics  have  to  date  been  limited  to 
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approximately  U  mm^.  It  has  been  suggested  that  the  cracking  in 
CdGeAs2  is  due  to  a  sphalerite  to  chalcopyrite  phase  transition  at 
g-jO°C!>3>39  We  have,  however,  some  doubt  about  the  phase  transition 
since  we  have  not  been  able  to  verify  the  sphalerite  phase  in  samples 
quenched  from  temperatures  above  630°C.  The  large  anisotropy  in  the 
thermal  expansion  (cf. Chapter  IV)  may  also  give  rise  to  cracking. 

There  is  a  definite  correlation  between  the  cracking  and  the  grain 
boundaries.  Large  single  crystal  regions  usually  have  few  cracks,  end 
it  is  hoped  that  the  cracking  will  reduce  with  improved  growth  techniques. 
Finally,  the  cracking  may  be  related  to  the  constraints  imposed  by  the 
walls  of  the  quartz  crucible. 

3.  Growth  Direction 

Wc  have  taken  Laue  photographs  of  several  boules  to  determine  the 
preferred  growth  direction.  In  a  few  boules  strain  caused  blurring  of 
the  Laue  spots.  The  orientation  of  large  single  crystal  regions  is 
usually  with  the  chalcopyrite  unit  cell  (111]  direction  within  10  of 
the  boule  axis.  This  agrees  with  the  literature  which  reports  the  [111] 
direction  as  the  fastest  growth  direction  for  the  chalcopyrites . 

The  impurity  segregation  coefficient  and  the  density  of  stoichio¬ 
metric  imperfections  may  depend  on  the  growth  direction.  Boule  number  31 
had  two  twins  with  a  large  disparity  in  the  optical  transmission.  The 
smaller  piece  was  transparent.  It  grew  in  a  direction  of  17  the 

c-axis  in  the  [100]  plane.  The  larger  piece,  however,  which  was 
completely  opaque,  grew  close  to  the  [221]  direction.  To  see  if  there 
was  any  significant  difference  in  the  impurity  concentration  between 
the  two  samples,  they  were  sent  out  for  impurity  analysis.  The  results. 
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however,  were  not  conclusive.  They  are  summarized  in  Table  X. 

4.  Optical  Transmission 

The  property  of  the  greatest  concern  is  the  optical  transmission. 
Since  the  bandgap  of  CdGeAs2  is  at  2. 3  pm,  visual  inspection  is  not 
possible.  Slices  1  nan  thick  are  therefore  cut  from  several  sections  of 
the  boule.  They  are  polished  and  the  optical  transmission  is  measured 
on  a  621  Perkin-Elmer  spectrophotometer.  Most  boules  have  a  very  non- 
uniform  transmission  with  the  best  transmission  usually  at  the  bottom. 
Uniform  transmission  exists  only  over  very  limited  regions.  The 
transmission  including  reflection  losses  through  one  of  our  best  samples 
of  p-type  CdGeAs^  is  shown  in  Fig.  13 .  More  often  the  best  transmission 
in  a  boule  is  between  20  to  30$  for  a  1  nm  thick  sample.  The  bandgap  at 
2.3  lira  and  the  two  phonon  absorption  at  18  pm  determine  the  transmission 
range.  Some  weak  three  phonon  absorption  exists  between  12  to  13  pm. 

There  is  no  free  carrie-  absorption  in  p-type  samples.  For  most  samples, 
however,  significant  absorption  occurs  between  the  bandgap  and  3  to 
6  pm.  This  absorption  is  not  seen  in  n-type  samples  which  in  return 
have  a  shorter  infrared  cutoff  wavelength  due  to  free  electron  absorption. 
Due  to  the  larger  bandgap  the  quality  of  CdGeP^  can  be  studied  using  an 
infrared  microscope  which  makes  evaluation  easier  than  for  CdGeAs^.  The 
boundary  between  transparent  and  opaque  regions  usually  follows  sharp 
lines  and  there  is  often  a  correlation  between  crack  lines  and  opaque 
regions  (cf.  Fig.  14).  The  cracking  is  concentrated  on  the  boundary 
between  single  crystal  regions.  The  optical  transmission  range  for 
n-type  CdGeP^  is  0.8  p  to  12.5  lira  as  shown,  in  Fig.  15.  There  is  some 


-  61  - 


(!N30d3d )  NOISSIINSNVdl 


-  62  - 


FIG.  15--Transmittance  of  a  single  crystal  of  p-type  CdGeAs 


FIG,  iU-“Transmlsslon  through  CdGeP0  photographed  through  an  IR  microscope. 
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Wavelength  ([in) 

FIG.  L5--Measured  transmission  range  for  n-type  CdGeP, 
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absorption  at  short  w.vrle  ngth.  close  to  the  bandgap  frequency.  Slmll.r 
absorption  is  seen  in  several  of  the  II-VI  semiconductor  compounds. 

5’  Measurement  of  the  Resistivity  and  the  Sian  of  the  MaW..,,  r . 

After  the  samples  have  been  tested  for  optical  transmission,  they 
are  etched  and  the  resistivity  is  measured  by  the  four  point  probe 
method.  In  addition,  a  hot  probe  is  used  to  determine  the  sign  of  the 
majority  carriers.  Growth  of  CdCeA.,,  from  stoichiometric  melt  gives 
P- type  material  except  for  a  region  close  to  the  c.pill.ry  where  there 
may  be  a  small  n-type  region.  Crystals  have  been  grown  from  melt  with 
0.5  to  1*  excess  arsenic,  and  one  boule  was  grown  with  2t  excess  ger¬ 
manium.  These  boules  were  also  p-type.  The  optical  transmission  end 
resistivity  of  boules  grown  from  arsenic  rich  melt  are  about  the  same 
as  for  stoichiometric  grown  boules.  Electron  probe  microanalysis 
revealed  araenlc  Inclusion,  near  the  top,  and  in  that  region  the 
material  was  „-type.  The  germanium  rich  boule  had  very  nonuniform 
transmission  due  to  germanium  precipitates. 

The  resistivity  correlates  with  the  optical  transmission.  For 
P-type  CdCeA.2  the  resistivity  at  room  temperature  1,  between  0.5(1  cm 
and  15!1  cm  with  the  largest  resistivity  resulting  in  the  best  optical 
transmission.  An  exception  was  boule  number  28  which  had  a  resistivity 
»f  about  100H  cm.  This  high  resistivity  was  probably  due  to  a  reduced 
mobility  because  of  strain  since  this  boule  was  quenched  from  i,80°C. 

The  presence  of  strain  was  confirmed  by  measuring  second  harmonic 
generation  in  a  wedged  sample. 
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The  resistivity  sometimes  varies  significantly  across  the  boule 
cross  section.  This  is  also  reflected  in  the  nonuniform  optical  trans¬ 
mission.  The  variation  in  resistivity  can  be  by  a  factor  of  five  or 
more.  For  example,  for  boule  number  31  which  consisted  of  two  single 
crystals,  one  •'ransparent  and  one  opaque,  the  resistivity  was  7.5  to 
15fi  cm  and  2.4  f l  cm  respectively. 

In  n-type  CdGeAs2  the  resistivity  is  usually  much  smaller  than 
1 n  cm  and  the  crystal  is  opaque. 

Only  a  few  boules  have  been  grown  of  CdGeP^  The  boules  were 
n-type  with  a  resistivity  of  10^  to  10^ a  cm. 

Ct  OTHER  TESTS 

1  -  Electron  Probe  Microanalysis 

We  have  made  extensive  microprobe  analysis  of  several  boules  to 
investigate  possible  deviations  from  stoichiometric  composition.  The 
absolute  accuracy  of  microprobe  measurements  are  between  2  and  5  weight 
percent.  For  absolute  calibration  we  use  Cd,  CdS,  Ge  ,  and  GaAs  as 
standards  and  determine  the  weight  fraction  Wu  of  the  elements  in  the 
unknown  (CdGeAs^)  using  the  expression 


W 


u 


(5-1) 


where  Wg  is  the  weight  fraction  of  the  element  in  the  standard  I 

’  s 

and  Iy  are  the  x-ray  ir tensities  from  the  standard  and  the  unknown, 
and  Fg  and  Fu  corrects  for  the  matrix  absorption  in  the  standard  and 
the  unknown. 
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We  have  determined  Fg/Fu  from  tables  in  Birks40  for  an  accelerating 
voltage  of  25  kV,  electrons  incident  at  62.5°  ,  and  a  takeoff  angle 
f  of  38-5°  •  Table  VI  gives  the  product  W  F  /F  .  Equation  (5.1) 
does  not  include  any  matrix  enhancement  factor  since  at  25  kV  there  is 
negligible  enhancement  in  the  x-ray  intensity  because  of  excitation  by 
the  characteristic  fluorescence  from  the  other  elements  in  the  matrix. 

Figure  16  shows  an  etched  boule  cross  section  of  CdGeAs  containing 
several  grains.  The  left  side  of  the  cross  section  was  transparent  and 
p-type  and  the  right  side  opaque  and  n-type.  Table  VII  gives  the  results 
of  the  microprobe  analysis.  No  variation  in  stoichiometry  over  the  cross 
section  or  irregularities  at  the  grain  boundaries  could  be  observed 
within  the  experimental  resolution. 

The  relative  accuracy  of  the  microprobe  analysis  can  be  as  good 
as  a  few  tenths  of  one  percent.  Table  VIII  lists  some  experimental 
results  for  different  boules  of  CdCeAs^  Within  the  experimental  error 
excess  As  does  not  perturb  the  stoichiometric  composition.  The  results 
agree  with  the  previous  phase  diagram  studies,  that  CdGeAs2  exists  only 
in  a  narrow  homogeneity  region. 

2.  Impurity  Analysis 

The  purification  of  compound  semiconductors  is  far  more  difficult 
than  for  the  element  semiconductors.  It  is  not  sufficient  to  remove 
foreign  atoms.  The  stoichiometric  composition  must  also  be  controlled. 

The  observed  carrier  concentration  is  probably  related  to  both  impurities 
and  stoichiometric  variations.  Since  there  are  approximately  Jo22  atoms/cm3 
a  carrier  concentration  of  lQ1^  cm  3  corresponds  to  only  1  ppm  of 
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TABLE  VI 


MAGNITUDE  OF  W  F  /F  FOR  MICROPROBE  ANALYSIS  OF  CdGeAs 
s  s'  u 


Standard 

Cd 

m 

Ge 

R2 

A3m 

Cd 

115.26 

CdS 

82.25 

Ge 

101.99 

GaAs 

1+9.57 

I  '  I  * 


TABLE  VIII 


MICROPROBE  ANALYSIS  OF  CdGeASg 


Boule  Number 

WCd 

w„ 

Ge 

CO 

** 

1 - 

1  Growth  Condition 

J _ _ 

Optical 

Transmission 

j'.6 

31.33 

1 

22.21 

46.46 

4 — - - — 

1 

' 0.5  Wt  $  excess  As 

Opaque 

;26A  (top  of  boule) 

31.33  | 

22.12 

46.55 

Stoichiometric 

Max  40$ 

26C 

31.34 1 

22.29 

46.36 

Stoichiometric 

i 

Max  15$ 

31 

31.64 1 

1 

22.03; 

46.34 

0.5  Wt  ^  excess  As 

1 

Max  10% 

33 

31.38 

I 

1 

22.36 

i 

46.26 

0.5  Wt  $  excess  As 

| 

Max  25# 

1 

CdGeAs 

2 

33.57 

21.68 

44.75  ! 

—  j 

1 

i 

— ! 

I 

j 

electrically  active  impurities  or  lattice  vacancies.  A  considerable 

amount  of  self  compensation  probably  occurs  in  CdGeAs  and  CdGeP 

2  2 ' 

Table  IX  shows  the  results  of  an  impurity  analysis  of  three 
slices  from  boule  number  33*  The  slices  were  from  different  sections 
of  the  boule  and  they  had  different  resistivity  and  optical  transmission. 
The  impurity  analysis  was  performed  by  the  Bell  and  Howell  Electronic 
Materials  Division  in  Pasadena  using  spark  source  mass  spectrometry. 

There  appears  to  be  very  little  correlation  between  the  impurity 
analysis  and  the  optical  transmission.  A  possible  explanation  is  that 
only  a  small  fraction  of  the  boule  cross  section  was  probed  and  homo¬ 
geneity  problems  may  have  obscured  the  results.  The  probed  area  was 
2 

less  than  one  mm  and  only  a  few  tenths  of  a  milligram  was  analyzed. 

The  usually  large  oxygen  and  carbon  concentrations  are  most  probably 
due  to  hydrocarbons  and  possibly  an  oxide  layer  on  the  surface. 

The  analysis  of  the  transparent  and  the  opaque  twin  in  Boule 
number  31  is  listed  in  Table  X.  Here  also  there  is  very  little  correl¬ 
ation  between  impurities  and  transparency.  Besides  oxygen  and  carbon, 
only  silicon  and  sulfur  impurities  are  present  at  a  significant  level. 

Most  of  the  crystal  growth  effort  has  centered  on  the  Bridgman 
method.  We  decided  to  investigate  this  method  since  it  is  a  relatively 
simple  method  which  usually  provides  good  results.  In  addition,  the 
method  had  already  been  used  by  others  to  grow  chalcopyrites.  The 
Bridgman  growth  program  is  now  near  completion.  We  have  been  unable 
to  solve  the  cracking  problem  and  the  optical  transmissica  \s  very 
nonuniform  throughout  the  boules.  Fortunately,  however,  the  necessary 
crystal  size  for  efficient  nonlinear  interactions  is  only  2  to  5  mm 
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and  from  some  sections  of  the  boules  we  have  obtained  crystals  of 
sufficient  quality  to  demonstrate  efficient  second  harmonic  generation 
of  a  C02  laser.  The  absorption,  however,  is  still  too  large  for  para¬ 
metric  oscillators.  We  have  tried  to  link  the  absorption  to  the 
impurities  or  stoichiometric  variations  using  mass  spectrometric  and 
electron  microprobe  analysis,  but  we  have  not  seen  any  significant 
correlation.  This  does  not  exclude  stoichiometric  variations,  however, 
since  the  microprobe  at  best  has  a  relative  accuracy  of  only  0.1*.  Much 
smaller  variations  in  the  stoichiometry  can  have  significant  influence 
on  the  absorption  and  carrier  concentric* .  With  reproducible  growth 
results  the  next  obvious  step  is  to  examine  how  compensation  and  also 
doping  affect  the  optical  transmission.  Alternative  growth  methods 
which  allow  seeding  and  avoid  the  constraints  set  by  the  crucible  wall 
should  also  be  explored.  With  a  large  anisotropic  thermal  expansion 
coefficient  c-axis  growth  is  probably  best  for  reducing  cracks.  It 
also  helps  to  grow  single  crystals.  In  the  Bridgman  boules  the  large 
strain  building  up  at  the  grain  boundaries  leads  to  cracking. 

Recently  we  have  tried  a  new  growth  technique  with  very  encouraging 
-esults.  By  growth  from  a  bismuth  solution  we  have  obtained  single 
crystals  with  almost  no  cracks.  The  optical  absorption  is  still  too 
high,  but  the  elimination  of  cracks  is  a  significant  step  in  the  right 


direction. 


TABLE  IX 


IMPURITY  CONCENTRATIONS  IN  CADMIUM  GERMANIUM  ARSENIDE 
(IN  PARTS  PER  MILLION  ATOMIC) 


Element^ 

Detection 

Limit 

Opaque 

Max  T  =  20$ 

Max  T  =  5$ 

Li 

0.007 

0.012 

O.O83 

0.017 

C 

0.03 

4.6-25 

3,200(75-10,000) 

19 

N 

0.03 

0.18 

0.51 

0.28 

0 

0.03 

76 

900 

3,600 

F 

0.07 

0.27 

0.1*8 

0.22 

Na 

O.Ol 

0.35 

8.3 

O.73 

Mg 

0.3 

N.D. 

0.1*9 

N.D. 

A1 

Si 

0.1 

1 

0.37 

6.7 

3-9 

N.D. 

0.80 

N.D. 

s 

0.03 

1-5 

1.1 

4.5 

K 

0.01 

0.052 

0.51  (6.6*) 

0.19 

Cd 

0.03 

O.O65 

0.25  (3.3*) 

0.071* 

(a) 


(b) 

* 

N.D. 


Waf  made  f°r  hydro8en’  Analyses  for  gold  are  not  given 
since  the  samples  were  sparked  against  high  purity  gold  countur- 
ectrudes.  Background  lines  of  the  matrix  interfere  with  the 

f°J  uC1’  and  Fe’  °ther  imPurities  not  listed  were  not 
detected  and  have  concentrations  less  than  0.3  ppma. 

etermined  for  3  x  10"7  coulomb  exposure. 

Seen  on  one  exposure  only. 

Not  detected. 
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TABLE  X 


IMPURITY  CONCENTRATIONS  IN  CADMIUM  GERMANIUM  ARSENIDE 
(IN  PARTS  PER  MILLION  ATOMIC) 


Boule  Number  31 


Element ^ 

Detection  Limit^*3) 

Transparent 

Sample 

Opaque 

Sample 

Li 

0.007 

O.O77 

0.021 

B 

C 

0.02 

0.0L6 

0.075 

0.03 

31 

110 

N 

0.03 

1.7 

2.0 

0 

0.03 

27 

120 

F 

Na 

Mg 

A1 

Si 

0.07 

0.01 

1 

0.1 

0.3 

0.1 

1.8 

N.D. 

1.8 

3.1 

0.1 

N.D. 

3.3 

16 

P 

S 

Cl 

K 

Ca 

0.03 

0.03 

300 

0.01 

0.03 

0.091* 

h.5 

N. D. 

1.5 

O. 78 

0.19 

1.8 

N.D. 

i*.7 

1.5 

V 

Cu 

Zn 

Ga 

Sn 

0.05 

0.07 

0.07 

0.1 

0.1 

0.20 

0.1*5 

0.2 

0.1 

O.65 

0.33 

0.93 

0.78 

0.29 

1.9 

I 

0.07 

0.39 

0.21* 

(a) 


(b) 

N.D. 


No  analysis  was  made  for  hydrogen.  Background  lines  of  the 

OthrVTr  With  the  analyS6S  m-ganese  and  iron. 

centraMPUri1 168  T  U8ted  W6re  n0t  detected  and  have  con¬ 
centrations  less  than  0.3  ppma. 

Determined  for  3  x  10"7  coulomb  exposure. 

Not  detected. 
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CHAPTER  VI 


LINEAR  AND  NONLINEAR  OPTICAL  PROPERTIES 

A .  THEORY 

In  this  chapter  we  briefly  review  some  nonlinear  optical 
theory  as  a  necessary  background  and  then  present  our  experimental 
results.  We  also  discuss  the  optical  absorption  mechanisms  in 
chalcopyrite  r/stals. 

1 .  Second  Order  Nonlinear  Interactions 
(a)  Effective  nonlinear  coefficient 
The  symmetry  restrictions  on  the  second  order  nonlinear  coefficient 
are  the  same  as  for  the  piezoelectric  tensor.  For  th-i  chalcopyrites 
which  have  42m  symmetry,  the  generated  nonlinear  polarizations  along 


the  principal  axes 

in  terms  of 

the  electric 

field  amplitudes  are 

6 

therefore 

P 

1  X 

=  2dl4EyEz 

P 

y 

-  2dl4EzEx 

(6.1) 

p 

z 

=  2d  /-E  E 

56  x  y 

• 

Since  ilh  =  dj6 

according  to 

the  Kleinman 

symmetry  condition,  the 

chalcopyrites  have  only  one  independent  nonlinear  coefficient.  The 
polarization  tensor,  which  is  the  same  as  that  for  KDP,  allows 
both  type  I  and  type  II  phasematching.  The  phasematching  or  conserva 
tion  of  momentum  conditions  for  the  type  I  and  type  II  phasematching 
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i 

f 


where  the  superscripts  e  and  o  denote  extraordinary  and  ordinary 

waves  and  p  ,  s  ,  and  i  refer  to  the  pump,  signal,  and  idler 

fields.  In  addition,  conservation  of  energy  requires  cu  =  cu  +  cu 

p  s  i 

For  type  I  phasematching  the  effective  nonlinear  coefficient  is 
dl4  sin  20  where  9  is  the  angle  between  the  direction  of  propagation 
and  the  c-axis.  The  signal  and  idler  are  extraordinary  waves  polarized 
in  the  (010)  plane  ind  the  perpendicular  pump  wave  is  ordinary.  Maxi¬ 
mum  nonlinear  interaction  occurs  at  0  =  45°  .  Only  type  II  phase¬ 

matching  allows  interaction  in  the  90°  direction.  The  effective 
nonlinear  coefficient  for  type  II  phasematching  is  d^  sin  0  with 
the  signal  polarized  as  an  extraordinary  wave  in  the(llO)  plane  and 
the  idler  polarized  parallel  to  the  pump  as  an  ordinary  wave.  Figure  17 
shows  how  the  fields  are  polarized  for  the  two  phasematching  conditions. 

For  second  harmonic  generation  (SHG)  the  Eqs.  (6.1)  reduce  to 


P2  =  di4  sin  29 

pl  =  2diu  sin  20  ele2  (6.4) 

for  type  I  phasematching  where  the  subscripts  1  and  2  refer  to 

the  fundamental  and  the  second  harmonic  waves .  The  fundamental  wave  is 
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A 


HpCUp  =  nj(0)cu8  +  n*(0)ui, 

Pp  *  2d,4sin20  EsEj 
Ps  »2d!4sin20E.Ep 
Pj  =  2d|4sin20EsEr 


['»] 


B 


FIG.  17 — Type  I  and  type  II  phasematching  in  a  positive  bire- 
fringent  crystal  of  T?2m  symmetry. 
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polarized  along  Eg  or  ^  in  Fig.  19(a)  and  the  second  harmonic  in 

the  E  direction.  For  type  II  phasematching  we  have 
P 

P2  =  dl4  sin  9 

P  =  2d^  sin  9  E^Eg  .  (6-5) 

Referring  to  Fig.  19(b)  the  fundamental  wave  is  now  polarized  in  the 
(E  E  )  plane  1+5°  to  the  E  direction  and  the  second  harmonic  along 

E  or  E, 

P  i 

(b )  Second  harmonic  conversion  efficiency 

Starting  from  Maxwell's  equations  we  can  determine  the  second  har¬ 
monic  conversion  efficiency.  This  is  a  standard  calculation,  and 
reference  4  gives  an  excellent  treatment  of  the  subject.  For  reference 
we  therefore  only  state  the  results.  Introducing  an  effective  non- 
linear  coefficient  d  anch  that  =  dB*  ,  «e  have  that  the  second 

4 

harmonic  power  is  given  by  the  equation 

P2  =  kLh  ,  (6.6) 


where  is  the  laser  power  at  the  fundamental  frequency  and  the 

constant  K  in  mks  units  is  given  by 


K 


5  2  2 

nln27T 


(6.7' 


The  dimensionless  quantity  h  in  Eq.  (6.6)  includes  the  effects  of 
double  refraction  and  focusing.  It  is  tabulated  in  reference  4  and 
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is  of  the  order  of  unity  or  smaller.  Further  f  is  the  crystal  length, 
nl  and  are  the  indices  of  refraction  at  the  fundamental  and 

second  harmonic  frequencies,  k^  =  n^cD^/c  is  the  wavevector,  and 
finally  T)q  =  377  f!  is  the  free  space  impedance.  In  order  to  compare 
non! inear  materials,  a  figure  of  merit  is  often  defined  given  by 


M 


1 


(6.8) 


This  definition,  however,  neglects  the  limitations  set  by  double 
refraction  for  phasematching  at  other  angles  than  90°.  To  correctly 
account  for  this,  we  introduce  two  other  figure  of  merits  which  apply 
to  the  cases  of  maximizing  the  total  SHG  power  and  the  SHG  power  in 
the  Gaussian  mode.  We  first  define  a  few  parameters.  The  fundamental 
wave  has  a  spot  size  w^  located  at  the  center  of  the  crystal  and 
the  confocal  parameter  is 

? 

27m.  w. 

b  -  •  (6.9) 

h 

The  parameters  E  and  B  describe  the  focusing  and  the  double  refrac¬ 
tion.  They  are  given  by 

.  i  1/2 

?  =  l/h  and  B  =  -  p(|kL)  ,  (6.10) 

where  p  is  the  double  refraction  angle  defined  in  Appendix  A.  We 
consider  first  the  case  of  total  SHG  power.  For  optimum  focusing 
(?m  =  1-392)  and  B  >  2  we  have  that  h  in  Eq.  (6.6)  reduces 
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to  the  approximate  expression 


hmm  =  °*71Vb  •  (6.H) 


From  Eqs .  (6.6),  (6.10),  and  (6.11)  it  follows  that  under  these  condi¬ 
tions  the  SHG  power  Improves  only  with  the  square  root  of  the  crystal 
length  and  the  appropriate  crystal  figure  of  merit  is  given  by  expression 


MSHG  tot  3/2 

nl  *i2p 

For  convenience  we  write  the  condition  B  >  2 
length  la  which  is  defined  by 


(6.12) 


in  terms  of  the  aperture 


i 

a 


(6.13) 


since  sometimes  this  parameter  is  used  instead  of  B  .  The  aperture 

length  should  be  longer  than  the  crystal  length  in  order  for  the  SHG 

power  not  to  be  limited  by  the  double  refraction.  With  B  =  (tt^2/2) 

{l/t  )?  1/^2  >2  we  obtain  l  <0.3  l 
'a  m  a 

For  many  applications  diffraction  limited  beams  are  of  importance. 
The  maximum  SHG  power  in  the  Gaussian  mode  is  obtained  by  substituting 
the  expression 


h  (B) 


b  (0) 

mm  ' 


mm 


1  +  (^B2/7r)hinm(o) 


(6. 1M 


into  Eq.  (6.6).  Here  h  (0)  is  equal  to  1.068.  Equation  (6.14)  holds 


to  within  10*  for  all  B  .  It  is  clear  from  Eqs .  (6.6)  and  (6.14) 
that  the  Gaussian  SHG  power  saturates  for  >  1  ,  and 

we  use  this  to  define  a  maximum  useful  crystal  length  imax  given  by 


9 

max 


(6.15) 


For  £  greather  than  imax  the  SHG  power  is 


and  in  this  case  the  figure  of  merit  is  given  by 

d2 

MSHC  Gaus  =  n2ng7 


(6.16) 


(6.17) 


We  later  use  the  figure  of  merit  expressions  to  compare  the  chalco- 
pyrltes  with  other  existing  Infrared  nonlinear  materials. 

For  weak  focusing  (E  <  1)  and  negligible  walk-off  (B  <  0.5),  »e 
have  h  ~  E  and  Eq.  (6.6)  reduces  to 


This  equation  Is  valid  for  small  conversion  efficiencies.  It  assumes 
that  the  phasematching  condition  Is  satisfied  and  neglects  absorption 
For  large  conversion  efficiencies  the  exact  expression  given  by 


(6.19) 
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must  be  used.  If  the  phasematching  condition  it;  not  satisfied, 
Eo.  (6-18)  modifies  to 


P, 


K 


*- 

P^  —  sine  [Hc.i/2) 

v„ 


(6.20) 


where  represents  the  momentum  mismatch .  In  the  above  equation;? 

P2  and  P1  are  the  powers  inside  the  crystal.  For  an  un-coated 
crystal  we  have  to  take  into  account  surface  reflections.  Using 
the  transformations  P, 


=  Ppl^iu/Cn  +  l)c]  and  p 


r  C. 

Pjl^l  -  1)  /kn^3  where  P£  _  _  and  P, 


1  out 
are  the  powers  out- 


out  ’  1  out 

side  the  crystal,  we  obtain  from  the  Eqs.  (6.7)  and  (6.20)  that 


P 


128  W<\  out 


2  out 


b  2 


(n2  +  l)  (n]L  +  1}  tiv2 


-a2//2  -a  i  -(a  /2  +  a  )/ 

(e  -  e  )2  +  he  1  sin2  £kl/2 

(&f  -i  (a2/2  -  at)2 


(6.21) 


In  this  last  equation  we  have  also  included  the  optical  absorption  at 
the  two  frequencies.  The  only  assumptions  in  Eq.  (6.21)  are  weak 
focusing  and  moderate  conversion  efficiency. 

(c)  Second  harmonic  angular  half-width 

Provided  that  the  absorption  is  not  too  large,  the  SHG  output 
power  from  a  crystal  follows  a  sine  curve  when  the  crystal  rotates 
through  the  phasematching  position  [cf.  Eqs.  (6.20)  and  (6 .21)]. 

Since  the  angular  half-width  depends  on  the  crystal  length,  we  can 
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use  the  observed  width  to  determine  an  effective  phasematching  length 
which  for  a  high  quality  crystal  is  equal  to  the  crystal  length. 

Variations  in  the  indices  of  refraction,  however,  may  make  it  impossible 
to  satisfy  the  phasematching  condition  over  the  total  crystal  length, 
leading  to  a  shorter  phasematching  length.  Since  sine2  80°  *  0.5  , 
it  follows  that  the  angular  half-width  at  the  half  power  point  is 

determined  by 

^  60  (6.2a) 

2  18c 

Writing  Ak(0 )  =  [3Ak(0)/Se]A0  ,  we  find  that 

kn  Bn® 

Ak{0)  =  - -  AS  (6-25) 

39 

for  type  I  and 

2n  Bn?  ..  .  . 

Ak(0)  =  - -  Ae  (6-24) 

xx  30 

for  type  XI  phasematching  In  a  positive  birefringent  crystal.  We 
determine  dn'/dS  by  differentiating  Eq.  (A.l)  in  Appendix  A  and  find 

— —  .  -  ne(0)  tan  p  .  (6.25) 

se 

By  combining  the  Eqs .  (6.22)  to  (6.25)  we  determine  the  internal  angular 
half-width  A0int  •  This  must  be  related  to  the  measured  external 
half-width  Aeext  .  When  the  laser  is  incident  on  the  crystal  at  an 
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angle  a  to  the  normal,  the  internal  angle  6  follows  from  the 
relation 


n1  sin  £  =  sin 


(6.26) 


If  the  crystal  phasematches  at  this  position,  we  immediately  obtain 


by  differentiation  that 


cos  f3 


(6.27) 


td.  . 

xnt 


cos  a 


For  close  to 


normal  incidence  this  equation  reduces  to  A0ext  - 


and  with  this  we  find 


1>ext  tan 


(6.28) 


n°  +  n*(9) 


II, ext  2.251  tan  p  2n®(0) 


< 6.29 ) 


for  the  t„o  phasematching  conditions.  The  half-widths  are  calculated 
in  radians  and  we  see  that  the  type  II  phaaematching  half-width  is 


twice  the  type  I  half-width 


In  our  experiments  we  have  mixed  5-5  m  with  10.6  nn  to  generate 
5.55  „m  using  type  I  phasematching.  For  this  experiment  the  angular 


half-width  is  given  by 


as5. 3+10.6  , 

I,  ext 


Xt  ln«(e)  +  Hg(e)] 


9i[ng(0)  tan  pg  +  1 /2  “^(9)  tan 


(6.50) 
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(d)  Parametric  oscillator  threshold  and  banawldth 


When  available,  CdGeAs^  is  very  useful  as  the  nonlinear  element 
in  an  infrared  parametric  oscillator  because  of  the  extensive  tuning 
range  and  large  nonlinear  coefficient.  In  this  section  we  list  the 
equations  which  allow  us  to  estimate  the  parametric  oscillator  threshold 
for  doubly  and  singly  resonant  oscillators.  In  a  doubly  resonant 
oscillator  (DRO)  there  is  low  loss  both  at  the  signa1  and  the  idler 

wave,  whereas  for  a  singly  resonant  oscillator  (SRO)  only  one  of  the 

waves  is  resonant.  The  single  pass  gain  of  a  degenerate  parametric 

amplifier  cjd  =  cjd.  -  cd  /2  is  equal  to  the  second  harmonic  conversion 
sip 

efficiency  Pg/P^  into  a  Gaussian  mode  when  the  fundamental  frequency 

is  at  cjd  /2  .  We  follow  reference  5  am*  define  cn^  =  CDp/2  , 

oj  =  a)  (1  +  8),  and  cjd  =  cjd  (1  -  &)  .  For  maximum  parametric  gain 
s  0  i  U 

both  the  pump,  signal,  and  idler  wave  should  have  the  same  confocal 

parameter  such  that  bp  =  bg  =  b±  .  This  leads  to  a  small  signal 

5 

gain  G  given  by 

G  =  (rf)2  =  KPpik0h(l  -  &2)2  ,  (6-51) 

where  K  is  defined  in  Eq.  (6.7)  with  =  cd^  .  Further  P  is 
the  pump  power  and  the  parameter  6  tells  how  close  the  oscillator 
is  to  degenerate  operation.  The  bar  expresses  that  only  the  coupling 
into  the  Gaussian  signal  and  idler  mode  should  be  included  in  the  gain 
expression  and  at  optimum  focusing  Eq.  (6.14)  gives  a  good  approximation 
for  h  •  It  is  clear  from  Eq.  (6.3I)  that  the  gain  decreases  when  the 
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osci] lator  tune  j  away  from  degeneracy, leading  to  an  increased  pump 
threshold.  With  single  pass  power  losses  ag  and  0^  at  the 
signal  and  idler  frequency  the  necessary  pump  power  to  reach  DRO 
threshold  is  given  by 

G  =  a  a.  ,  (6.32) 

si 


and  for  SRO  operation  with  no  idler  feedback  (a^  =  l)  we  have 

G  =  2a  .  (6.33) 

s 


These  equations  are  strictly  valid  only  for  cw  operation.  For  a 
Q-switched  pump  laser  the  threshold  must  be  reached  within  a  finite 
number  of  passes,  and  this  requires  some  excess  gain.  Assuming 
a  square  pump  pulse  of  length  t  and  a  cavity  transit  time  t^ 
the  maximum  number  of  available  passes  is  n  =  t^Ag  •  Under 
these  conditions  it  is  easy  to  show  that  the  DRO  cw  threshold  condition 

41 

modifies  to 


/r  =  ■  a  * 


(6.310 


when  the  signal  and  idler  wave  has  the  same  single  pass  power  loss  a  . 
This  equation  gives  the  necessary  gain  for  the  oscillator  to  build  up 
from  the  parametric  fluorescence  noise  power  P^  at  the  signal  or 
idler  wavelength  to  a  power  Pn  after  n  passes.  We  require  some 
pump  depletion  and  take  Pn  to  be  10$  of  the  pump  power.  Because  of 
the  logarithmic  dependence  the  exact  value  of  Pn/pQ  not  critical. 
As  an  example,  for  p^/p^  «  10  L  ,  t  =  200  nsec,  and  t^  =  0.2  nsec 
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which  corresponds  to  an  effective  cavity  length  of  6  cm,  the  equivalent 
loss  due  to  buildup  is  1  .bio.  Similarly  for  pulsed  SRO  operation  we 

k-2 

modify  Eq.  (6.33)*  Resonating  the  signal  wave  leads  to 


1 

lg( 1  +  G  -  2a g)  =  -  lg 

n 


(6-35) 


assuming  no  loss  in  the  crystal  at  the  idler  wavelength  and  also  that 


G  -  2  «  1  .  The  eouation  reduces  to 

x 

s 


G  =  2a  + 
s 


?  -  § 


(6.36) 


Using  the  same  numerical  example  as  before  the  necessary  excess  gain 
is  2. bio. 

2 

The  frequency  output  of  a  parametric  oscillator  follows  a  sine 

(Aki/2 )  dependence  where  Ak  =  k  -  k  -  k  is  the  momentum  mis- 

p  S  1 

match.  Assuming  a  single  frequency  pump,  a  change  6u)  (=  -Sen  ) 

s  1 

in  the  signal  frequency  away  from  the  phasematching  solution  leads 

1+1 

to  a  momentum  mismatch  or 


where 

b 


Ak  =  b5<ug  +  |  bj^  ,  (6-37) 


(6.38) 
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and 


3b 


Qu: 


2nc 


.2  .2 

an  ,  5  ni 


(6.39) 


We  define  the  bandwidth  by  setting  |Aki/2|  =  7T  .  Except  for  operation 

close  to  the  degeneracy  point  for  type  I  phasematching  or  near  the 
turning  points  of  the  angular  tuning  curves  for  type  II  phasematching 
we  have  |  krtoJVo2 1<1  and  by  solving  Eq.  (6. 37)  for  6mg  we 
determine  the  full  bandwidth  Am  at  the  .base  line.  We  find 


For  the  special  case  that  |  Imb^/Lb1^  |  >  1  ,  the  bandwidth  is  given 


by 


(6.41) 


2 .  Third  Harmonic  Generation 

For  an  electron  concentration  of  5  X  10^  cm  ^  the  III-V  compounds 

43 

InSb  and  InAs  have  very  large  third  order  nonlinear  coefficients 
with  the  main  contribution  to  the  nonlinearity  arising  from  the  non- 
parabolicity  of  the  conduction  band.  The  magnitude  is  inversely 
proportional  to  the  bandgap  frequency  and  the  effective  mass  squared 
and  proportional  to  the  carrier  concentration.  This  makes  CdGeASg 
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interesting  because  in  addition  to  the  .11  bandgap  frequency  and 
the  small  effective  mass  it  has  sufficient  birefringence  for  phase¬ 
matching  leading  to  the  possibility  of  reasonably  efficient  THG.  A 
determination  of  the  various  components  of  the  susceptibility  tensor 
may  also  provide  an  interesting  check  on  the  bandstructure  calculation 
neat  k  =  0  .  In  this  section  ve  discuss  phasematching  and  derive 

expressions  for  the  effective  nonlinear  coefficient,  conversion  effi¬ 
ciency,  and  angular  bandwidth  for  third  harmonic  generation  (THG). 

(a)  thp.  phase-match1' np  and  effective  nonlinear  coefficient 

There  are  three  possible  ways  to  phasematch  THG.  For  a  crystal 
with  positive  birefringence  the  phase, .etching  conditions  can  be  written 


as 


>5 


n°  =  nA9) 

3  1 

n°  =  i  [2n®(e)  +  n°) 
n°  =  i  ln*(0)  +  2n°] 


(6.1+2) 


and  we  refer  to  them  as  respectively  type  I,  It,  and  III  phasematching. 
Type  III  requires  the  largest  birefringence. 

We  use  cljk,  for  the  third  order  nonlinear  coefficient.  The 


definition  follows  from  the  equation 


?i  =  Cijk£EjEkE£ 


(6.1+3) 


with  summation  over  repeated  indices  and  where  as  usual  P  and 
stand  for  the  generated  nonlinear  polarization  and  the  electric  field 
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amplitudes.  The  permutation  symmetry  between  the  last  three  indices 

suggests  introduction  of  a  contracted  notation.  The  matrix  cljkf 

.  ..  „  hn<?  onlv  7>0.  The  definition 

has  81  elements  while  the  compact  form  cim  has  only 

of  m  follows  from  Table  XI.45  For  the  42m  chalcopyrite  symmetry 
the  tensor  cim  is  given  by 

0  cl6  0  cl8  °  ° 

o  o  o  0  ci8  0 

.  0  C  o  o  0 

35  55 

(6.44) 


c . 
im 


'11 

0 

0 


0 

'11 

o  . 


0 

0 

c33 


:l6 


The  tensor  has  5  independent  elements  or  4  when  we  use  the  Kleinman 


symmetry  condition  which  gives  c35  =  C16 


Eq.  (6.43)  reduces  to 


Pi  ~  cimSm  * 


where  the  vector  is  given  by 


111 


45 


Sm  " 


222 


333 


+  £ir. 


233  323  532 
^223  +  ^232  +  ^322 
^133  +  Sl5  +  ^31 
^113  +  ^131  +  ^311 
^122  +  ^212  +  ^221 
*112  +  *121  +  *211 


In  contracted  form  the 


(6.45) 


'123 


+  *132  +  *213  +  *231  +  *312  +  *321 


(6.46) 


TABLE 


XI 


CONTRACTED  NOTATION  FOR  THE  THIRD  ORDER  NONLINEAR  COEFFICIENT 


i  jk  1*1 

■ - 1 - 

222 

333 

233  223 

1 

1 

133  113 

122 

jgm 

m  J  1 

1 

1 

2 

3 

i 

1 

^  5 

i 

1 

■ 

6  ;  7  j 

■ 

l  i 

8 

mi 

anG  ^ijk  -  (^2 ^1  ^  ‘  use  •  (6.45)  together  with 

the  Eqs .  (6.44)  and  (6.46)  to  determine  the  effective  nonlinear 


coefficient  for  the  three  phasematching  conditions.  Let  9  be  the 
angle  between  the  wavevector  and  the  z  axis,  <p  the  angle 

between  the  projection  of  ^  in  the  (x,y)  plane  and  the  x  axis, 
and  a  the  angle  between  the  electric  field  at  the  fundamental 

frequency  and  the  normal  to  ^  in  the  (k^z)  plane  (cf.  Fig.  18). 
The  effective  nonlinear  coefficient  depends  on  the  angles  and  we 
find 


I: 

P3 

■  -Uhl' 

II: 

P3 

=  [2  (Cll 

2  1 

+  c^  cos  e 

III: 

P3 

•  i  (cn  ' 3c 

32a  eJ 


(6.47) 


for  the  three  phasematching  conditions.  In  Eqs.  (6.47)  9  is  fixed 

by  the  phasematching  conditions  in  Eqs.  (6.42).  The  two  other  angles, 
however,  can  be  chosen  to  maximize  the  nonlinear  interaction.  The 
optimum  values  for  cp  and  a  are  listed  in  Table  XII  together  with 
the  maximized  effective  nonlinear  coefficient  c  as  defined  by 
^  =  c^]_  •  Considering  only  the  nonparabolic  band  contribution  to 

the  nonlinear  coefficient  we  have  for  spherical  bands  that 
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**y 


FIG.  l8--Definition  of  the  angles  6,  qr,  and  a  for  third 
order  optical  interactions. 
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i 

! 


-c  =  c  /-  =  c,o  .  With  this  assumption  the  effective  nonlinear 

5  11  16  lo 

coefficients  in  Table  XII  reduce  to  Cj  =  cnI  =  0  and  cn  = 

Since  the  chalcopyrites  have  fairly  isotropic  conduction  bands,  the 
most  efficient  third  order  interaction  in  n-type  material  probably 
occurs  for  type  II  phasematching. 

(b)  THG  conversion  efficiency  and  angular  bandwidth 

For  weak  focusing  such  that  the  crystal  length  l  is  smaller  than 
the  confocal  parameter  b^  ,  the  third  harmonic  conversion  efficiency 
is  given  by 

'i  _  >  > . i&  f  ,c.m 


,  {CM) 


2 

where  c  is  the  effective  nonlinear  coefficient  and  Aj  =  -nWj/2  = 
^b^/4n^  is  the  beam  area.  Equation  (6.48)  does  not  include  reflection 


losses . 


The  sine  curve  determines  the  angular  half-width.  Similar  to  the 
treatment  of  SHG  we  find  that  the  angular  half-widths  at  the  half 
power  point  are  given  by 


!,ext  tan  p 


II,ext  4.5 /  tan  p  n®(e) 


III, ext 


2.25 1  tan  p  n®(e)) 


(6.49) 


3 .  Determination  of  Indices  of  Refraction 


In  order  to  determine  accurately  the  phasematching  conditions  for 
a  material,  the  dispersion  and  the  birefringence  must  be  known  to  better 
than  one  part  in  the  third  place.  We  measure  the  indices  of  refraction 
using  a  prism  set  of  minimum  deviation.  This  is  the  best  method 
when  sufficiently  large  single  crystals  are  available. 

Figure  19  shows  a  diagram  of  the  index  of  refraction  measurement 
apparatus.  The  light  source  is  either  a  laser  or  a  globar  for  measure¬ 
ments  at  wavelengths  longer  than  one  micron.  The  detector  is  FbS  out 
to  three  microns,  and  then  a  thermocouple  at  longer  wavelengths.  ine 
crystal  is  mounted  on  a  Gurley  Unisec  table  which  measures  angles  to 
one  second  of  arc.  In  practice,  diffraction  effects  due  to  the  finite 
prism  size  limit  the  accuracy  of  the  measured  angles.  With  the  col¬ 
limated  light  filling  a  prism  of  length  L  ,  Fraunhofer  diffraction 
limits  the  angular  width  of  the  focused  beam  at  the  detector  to 

Xf 

&P  -  2  —  ,  (6.50) 

LR 

where  X  is  the  wavelength,  R  is  the  distance  from  the  center  of 
the  table  to  the  detecto. ,  and  f  is  the  focal  length  of  the  focusing 
optics.  For  our  system  we  have  R  -  23  cm  and  f  =  15  cm.  Assuming 
the  detector  can  be  set  vo  the  maximum  within  five  percent  of  the  full 
angular  width,  we  have  for  a  0.5  cm  prism  and  a  wavelength  of  5  (Im  that 
the  diffraction  limits  the  accuracv  of  the  measured  angles  to  approxi¬ 
mately  14  seconds  of  arc. 
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4  > 


We  calculate  the  indices  of  refraction  from  the  equation 


n 


air 


sin  (a/2) 


(6.51) 


where  &  is  the  minimum  deviation  angle  and  a  is  the  apex  angle. 

By  differentiating  Eq.  (6.5I)  we  obtain 

sin  (6/2) 

A - - - 62  • 

2  sin  (a/2)  sin((a  +  &)/2) 

(6.52) 

For  our  CdGeAs2  and  CdGeP2  prisms  we  have  &  *  35°  «nd  a  * 

Suljtitution  into  Eq.  (6.52)  yields 

—  *  1.1  A&  -  3-3  &  >  (6*55) 

n 

where  A5  and  £&  are  in  radians,  and  the  uncertainty  in  the 
measured  index  values  is 

An  =  n  yj^. 1  A&)2  +  (3-5  tnf  -  (6-5*0 

Taking  A&  =  ±  1'  and  to  •  ±  V  the  uncertainty  is  An  =  ±  0.004  . 
The  absolute  accuracy  of  the  measured  indices  is  probably  not  better 
than  this  since  we  worked  with  very  small  prisms.  More  important, 
however,  is  the  relative  accuracy  which  depends  only  on  A5  .  We 
expect  this  to  be  better  than  one  part  in  the  third  place. 

The  sample  temperature  must  be  kept  constant  during  the  measure- 
M„t,.  For  cdCe*s2  dn/dc  ~  5  x  10'k  V1  .t  5-59  o»  •  This 
that  a  temperature  rise  of  2°C  changes  the  indices  by  as  much  as  one 
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part  in  the  third  place.  The  growth  method  and  carrier  concentration 
must  also  be  specified.  In  semiconductors  carrier  concentration 
variations  may  change  the  indices  of  refraction  significantly.1*7 
If  large  enough  single  crystals  are  not  available  for  prism 
fabrication,  other  methods  can  be  used  to  determine  the  birefringence 
of  a  crystal .  One  way  is  to  use  a  laser  to  probe  a  wedged  platelet 
between  crossed  polarizers.  The  laser  beam  is  normal  to  the  crystal 
and  polarized  1*5°  to  the  "effective"  c-axis  in  the  platelet  plane. 

By  translating  the  wedge  a  distance  x  across  the  beam,  the  power 
P  at  the  detector  is  given  by 


where  r,  -  T0  =  (2ff/x)[ne(fl)  ■  nj  x  Un  Q  is  the  phase  difference 
between  the  extraordinary  and  the  ordinary  wave.  Here  0  as  usual 
is  the  angle  between  the  c-axis  and  the  propagation  direction  and  a 
is  the  wedge  angle.  If  the  distance  between  two  power  minima  is  x 

o  f 

then  the  equation 


x  tan  a 
o 

determines  the  birefringence. 

For  very  small  samples  not  larger  than  1  mm,  the  birefringence 
can  be  measured  by  rotating  a  platelet.  This  method,  however,  is  not 
so  reliable.  With  crossed  polarizers  and  the  laser  polarized  at  hf 
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to  the  "effective"  c-axis,  Eq.  (6.55)  still  holds,  and  the  phase 
difference  is  given  by 

271 7 

re  -  r0  =  -  [n£(e)cos  cpe  -  nQ  coscpo)  ,  (6.57) 


where  t  is  the  plate  thickness  and  Fig.  20  defines  the  angles. 
Equation  (6.57)  can  be  written  as 


To  obtain  an  approximate  expression  we  use  Eq.  (A  .2)  in  Appendix  A 

2  2 

and  take  nQ  »  sin  ®  .  This  gives 


27! t  p  /  1  sin^cp 

re  -  ro  -  —  4n  sin  0  (l  ♦ - — 

1  \  2  no 

where  An  is  the  birefringence.  When  the  rotation  is  around  the  axis 
normal,  to  the  "effective"  c-axis,  9  is  related  to  rhe  external 
angle  q>  by  9  =  7r/2  +  9  -  q  where  q  is  the  angle  between  the 
c-axis  and  the  platelet  plane.  By  measuring  the  external  angle 
difference  between  two  power  minima,  we  use  Eq.  (6.55)  together  with 
Eq.  (6.58)  or  (6.59)  t0  calculate  the  birefringence. 

The  measured  index  data  can  be  fit  to  a  classical  Sellmeier 
equation  of  the  form 


(6.59) 


2 

n 


A  + 


(D/x)2  1  -  (E/x)2 


(6.60) 
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which  holds  in  the  transparent  region  of  the  crystal  not  too  close  to 
the  absorption  bands  .  In  this  equation  E  is  approximately  the 
reststrahl  wavelength  and  D  corresponds  to  an  average  electronic 
bandgap  wavelength.  From  Eq.  (6.60)  we  estimate  the  DC  dielectric 
constant  by  letting  X  -*  00  •  This  gives 

e__  =  A  +  B  +  C  .  (6.61) 


1+.  Measurement  Technique  for  the  Optical  Nonlinear  Coefficient 


There  are  several  ways  to  measure  the  optical  nonlinear  coefficients . 
We  have  used  the  wedged  sample  technique**  which  is  most  convenient  for 
relative  measurements.  Translation  of  the  wedged  nonlinear  crystal 
across  the  l.iser  beam  results  in  oscillations  in  the  SHG  output  power. 
Taking  s  as  integer,  maxima  occur  every  time  the  effective  crystal 
thickness  is  equal  to  2s  +  1  times  the  cohe.ence  length  which 

for  SHG  is  defined  as 


7T  Xj^ 

Ak  U(n2  -  nx) 


(6.62) 


1+9 

The  wedge  method  is  very  similar  to  the  Maker's  fringe  technique, 
but  it  has  some  advantages.  The  Maker's  fringe  technique  is  only 
useful  when  the  crystal  length  t  >  (2n)  / c  ,  and  the  method 
therefore  requires  large  samples  of  materials  with  large  indices 
of  refraction  and  a  large  coherence  length  when  only  small  samples 
are  available.  The  wedge  technique  does  not  have  this  limitation. 
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Thin  samples  may  be  used  and  the  absorption  can  therefore  most  often 
be  neglected.  The  wedge  technique  also  avoids  problems  of  multiple 
reflections  and  resonance.  The  fraction  of  the  SHG  power  that  suffers 
multiple  reflections  is  transmitted  in  different  directions  and  can 
be  stopped  by  proper  aperturing.  By  mounting  the  wedge  such  that  the 
incident  laser  beam  is  normal  to  the  crystal  surface,  there  is  no 
mode  distortion  due  to  refraction  as  it  will  be  with  the  Maker's 
fringe  technique. 

Equation  ^B.5)  in  Appendix  B  gives  the  SHG  output  power  generated 
in  a  wedged  sample.  Neglecting  the  absorption  the  equation  reduces 
to 


128  r?Q  u)2  d2  P? 


lout 


2out 


(°2  +  1)  (nt  +  1> 


|4  7TW2 


cos(Ak/Q)  exp 


fWj  Ak  tan  a' 


(Ak)2/2 


(6.63)  J 


which  agrees  with  reference  50.  We  determine  the  coherence  length 
*c  *  "n/Ak  by  translating  the  crystal  a  distance  y^  between  two 
SHG  minima.  This  gives 

lc  =  2  y0  tan  Q  ‘  (6.6k) 

The  expression  inside  the  parenthesis  in  Eq.  (6.63)  reduces  to  the 


standard 


form  when 


tin. 

k 


l 


This  condition  is  never  difficult  tj  satisfy  and  neglecting  the.  expo¬ 
nential  we  have 


2out 


128  TjJ  GO?  d^ 

‘Oil  lout 

/  .  v  2  /  -  %  2 

(n2  +  1)  (nj  +  1)  TTWj^ 


2  ,  .  7T  *0 

-  /  sin  —  — 

7T  C  2  i 

c 


(6.65) 


We  measure  the  relative  effective  nonlinear  coefficient  by  using  a 
reference  sample  and  comparing  the  coherence  lengths  and  the  aaximum 
SHG  output  powers  for  a  given  input  power. 

5.  Absorption  Mechanisms 

A  maximum  optical  power  density  in  the  material  eventually  limits 
the  efficiency  of  nonlinear  optical  processes.  Different  mechanisms 
may  be  responsible  such  as  nonlinear  absorption,  broken  phasematching 
condition  due  to  induced  index  of  refraction  inhomogeneities, 
or  permanent  crystal  damage.  Both  the  broken  phasematching  condition 
and  the  crystal  damage  are  often  caused  by  the  local  heating  due  to 
the  optical  absorption,  and  it  is  therefore  important  to  reduce  the 
absorption  to  a  minimum.  Some  of  the  most  common  absorption  mechanisms 
are: 

1)  Bandgap  absorption, 

2)  Two-phonon  absorption, 
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5)  Three-phonon  absorption, 

4)  Free  electron  scattering;  in  n-type  semiconductors, 

5)  Free  hole  scattering  ) 

V  in  p-type  semiconductors, 

6)  Intra  band  absorption) 

7)  Impurity  and  defect  absorption, 

8)  Free  carrier  absorption  generated  by  impurity  or  defect  absorption, 

9)  Two -photon  absorption, 

10)  Free  carrier  absorption  generated  by  two-photon  absorption. 


As  discussed  earlier  the  direct  or  indirect  bandgap  absorption 
determine  the  short  wavelength  cutoff  and  two-phonon  absorption  limits 
the  transmission  at  long  wavelengths.  In  CdGeASg  the  two-phonon  absorp 
tion  is  approximately  20  cm"1.  Some  weak  three-phonon  absorption 
exists  at  wavelengths  roughly  two  thirds  of  the  two-phonon  cutoff  and 
in  CdGeAs2  the  magnitude  is  approximately  0.3  cm  at  12  pm. 

The  free  carriers  in  a  semiconductor  affect  both  the  indices  of 
refraction  and  the  absorption.  According  to  classical  Drude-Zener 
theory  the  free  carrier  contribution  to  the  relative  dielectric 


constant  e  at  frequency  co  is  given  by 


51 


=  € 


1  -  i 


CO  COpT 


CD  1  +  iCOT 


(6.66) 


where  e  is  the  optical  relative  dielectric  constant  without  the 

CO 

carrier  contribution,  co  the  plasma  frequency,  and  t  the  carrier 

P  * 

collision  time.  For  a  carrier  density  N  and  an  effective  mass  m 
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the  plasma  frequency  is  given  by 


(6.67) 

In  the  visible  and  near  infrared  where  co  »  oip  and  1/t  ,  it  follows 
from  Eq.  (6.66)  that  the  indices  of  refraction  n  and  the  free  carrier 
absorption  constant  a  can  be  written  as 


and 


n 


!(u>p/a>)2 


(6.68) 


a 


2  2 
n  os 

oo  p 


C  I  (11 


(6.69) 


In  crystals  with  an  anisotropic  effective  mass  tensor,  the  presence  of 
carriers  changes  the  birefringence.  Assuming  a  single  band  extrema 
at  k  =  0  ,  the  free  carrier  induced  change  in  the  birefringence  is 
found  from  Eqs .  (6.67)  and  (6.68)  to  be 

1  1 

(6.70) 

We  refer  to  this  expression  in  Section  A. 6  of  this  chapter  and  discuss 
how  to  use  this  effect  in  CdGeAs^  to  construct  a  modulator. 

According  to  Eq.  (6.69),  the  free  carrier  absorption  increases 
with  the  square  of  the  wavelength.  This  result  assumes  an  energy 
independent  collision  time  which  often  is  not  the  case.  More  generally, 


Ne 


nii  -  n  1  =  mi 


“1 


-  n  | 

,°0  _[_,<*> 


2e0  “ 


we  can  write  a  *  x"  where  the  exponent  depends  on  which  scattering 
mechanism  limits  the  collision  time.  Scattering  by  acoustical  phonons 
gives  1.5, 52  and  scattering  by  optical  phonons  gives  2.5.  For 

5I1 

ionized  impurity  scattering,  the  exponent  is  between  3.O  and  3.5, 

depending  on  which  approximation  is  used.  For  CdGeAs^  we  have  measured 

a  3.5  dependence  for  n-type  material  (cf.  CH.  VI.B.5).  The  absorption 

is  therefore  determined  by  ionized  impurity  scattering  and  the 

. ,  54 

absorption  cross  section  follows  from  the  equation 

16  Vi"  7T2  /  Ze2  \2  e2  N±  ft2 

°abs  “  *^el  3  %£0  n»  c^)3//2 

(6.71) 

where  N  is  the  density  of  ionized  impurities.  The  derivation  of 
Eq.  (6.71)  assumes  an  isotropic  effective  mass,  but  the  expression 
still  can  be  used  to  estimate  the  ionized  impurity  density  in  CdGeAs^ 
since  the  conduction  band  is  not  too  anisotropic. 

In  p-type  CdGeAs2  we  see  no  evidence  of  free  hole  absorption. 

The  absorption  has  a  more  complicated  wavelength  dependence  due  to 
intra  band  transitions  between  the  split  valence  bands .  This  absorp- 

CC 

tion  also  occurs  in  p-type  germanium^'  and  the  III-V  compounds. 

We  refer  to  the  results  of  Chapter  III  to  determine  the  magnitude 
and  spectral  dependence  of  the  intra  band  absorption  in  p-type  CdGeASg. 
The  absorption  constant  at  the  frequency  u>  for  transitions  between 
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the  v^  band  the  the  v£  or  v^  band  is  * 

E. /kT  , 

M|2  e  ^  B(E1  -  Ei  -  *to)d^k 

(6.72) 


—  1-  (i .  e'to/kT; 
.  2  n  cur  M 

47re0m  c  co  N 
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Vhere  1  =  2  or  3  refers  Co  the  two  lower  valence  bands,  m  Is  the 

free  electron  mass,  P  Is  the  hole  density,  and  Is  defined  In 

E,.  (b.9).  The  integration  Is  over  the  crystal  momentum  with 

i,  _  v  -  k  .  Further  the  matrix  element  M  Is  given  by 
K  “  K1  "  i 


|M|J 


A  “I  >A  t.2  i  /  .A 


)l2  +  l<<  J* 


pk®  )l‘ 


(6.75) 


and  the  energy  «,  of  the  v,  band  Is  measured  relative  to  the  band 
extrema  such  that 


.  ,2 


E, 


1T 

2 


Lmi,i  -i.i 


(6.74) 


we  evaluate  the  Integral  In  Eq.  (6.72)  in  Appendix  C.  Equation  (5.2b) 
gives  the  matrix  elements  for  the  Intra  band  transitions.  By  combining 
Eq.  (3.2b)  with  Eqs .  (C.3),  (C.12),  (C-U),  »”d  (6.72)  «e  calculate 

,  .  2  _  2  _  2  _  i  .  The  expression 

the  absorption  constant  letting  ex  y  "  z  5 

for  the  absorption  constant  depend,  on  only  four  parameters:  the  b.nd- 

g.p  energy,  the  spln-orblt  coupling,  the  crystal  field  splitting,  and 

the  matrix  element  P  .  With  the  values  In  Table  II  anu  Bq.  (5-20). 

a.  Crow  rHfpAs  at  a  carrier  concentration  of 
the  absorption  constants  for  CdGeASg  a 

2  x  lO1^  cra”^  can  be  written 

_L  fl  -  exp(-WkT)]  (0.616  -  0.956  eAT)  exp(0.088  e/Vl) 

a2,l  “  «d  1 

(6.75) 

_  i.  {1  .  exp(-itoAT)Hl.85  -  2.62  eA*)  exp(0.171  eA*) 

*5,1  “  «o  1 

(6.76) 
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For 


.here  €  -  AE  -  ito  <X)  and  AE  Is  the  band  splitting  energy. 
the  Vi  .  v2  transition  .e  have  AE  =  0.20  eV  and  for  the  v,  - 
transition  AE  =  0.52  aV  .  Similarly,  .hen  t>0  »e  obtain 

n  -  exp ( -fiu)/kT ) ]  jlooljr  exp(-eAT) 

a2,l  ~  »» 

+  a,  ,  [1  -  erf  (/Till  eAU)  (f"77) 

Upl¬ 


and 


a 


,  _  n  -  exp(-WW))  /30o77r«p(-1.05  GAT)1  » 

3,1 

(6.78) 


.  a,  ,  |1  -  art  u  >■«  ^T)) 
Jr 


.here  a,  ,  and  a}  ,  are  given  by  the  same  expressions  as  the  .  < 
ehsorpti;:  constants  in  E,s.  (6.75)  and  (6.76).  Figure  21  sho^  the 
wavelength  dependence  of  the  c-ined  absorption  constants  for  the  two 
mtr.  band  transitions  at  three  different  temperatures.  For  a  carrier 
concentration  of  2  x  10*  -  ’  there  is  significant  absorption  .ithln 
the  normally  transparent  frequency  range.  The  long  wavelength  cutoff 
0f  o.20  eV  correspond,  to  the  calculated  splitting  energy  between  the 
V  and  v2  band.  In  the  experimental  section  of  this  chapter  we 
compare  measured  and  calculated  absorption.  The  observed  absorption 
cut-off  provides  an  experimental  value  for  band  splittihg  energy. 

Absorption  also  occurs  due  to  impurities  and  defect,  and  the 
structure  of  the  Intraband  absorption  is  therefore  often  obscured. 
Defects  such  as  lattice  vacancies  have  an  effect  similar  to  Impurities. 
The,  give  rise  to  acceptor  and  donor  levels  and  optical  absorption 


-  no  - 


Wavelength  (|iml 


m 


-  Ill  - 


Photon  energy  (eV) 

FIG.  21 _ Calculated  total  absorption  coefficient  for  transitors  between  the  split 

valence  bands  in  p-type  CdGeASg  at  three  different  temperatures  but  assuming 
a  constant  carrier  concentration. 


is  possible  between  these  levels  and  the  valence  or  conduction  band. 

If  the  imourity  levels  are  ionized  by  photon  absorption,  the  resulting 
carriers  lead  to  additional  free  carrier  absorption.  This  is  a  non¬ 
linear  absorption  mechanism  which  is  important  at  high  optical  power 
densities.  In  some  materials  nonlinear  absorption  leads  to  satura¬ 
tion  in  the  nonlinear  conversion  efficiency  before  the  optical  damage 
threshold  is  reached. 

Two-phctor  absorption  between  the  valence  and  the  conduction  band 
is  another  important  nonlinear  absorption  mechanism.  For  a  parametric 
oscillator,  the  pump  frequency  should  be  chosen  to  be  less  than  half 
Che  bandgap  frequency  to  avoid  two-photon  absorption.  The  measured 
two-photon  absorption  constant  in  GaAs  and  InP  is  approximately 
0.1  cm_1/taW  cm"‘  at  1 .06  pm.*57  The  free  carriers  generated  by  two- 
photon  absorption  lead  to  rn  absorption  constant  propo-tional  to  the 
square  of  the  laser  power  density.  In  InP  the  two-photon  generated 
free  carrier  absorption  equals  tha  normal  two-photon  absorption  at 
2  MW/cm2.  This  mechanism  is  believed  to  limit  the  10.6  iim  doubling 
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efficiency  in  tellurium  to,  at  best,  a  few  percent.  ’ 

6.  Modulation  bv  Free  Carrier  Induced  Birefringence 

Because  of  the  anisotropic  effective  masses  in  CdGeASg,  it  is 
possible  to  modulal  the  birefringence  by  injecting  free  carriers. 

The  maximum  modulation  frequency  is  eet  by  the  free  carrier  recombi¬ 
nation  time.  Increasing  the  carrier  concentration  by  AN  changes 
the  birefringence  by  6 (An)  which,  according  to  Eq.  (6.70),  may  be 
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written  as 


1 


{(An) 


assuming  the  injection  generates  equal  numbers  of  holes  and  electrons. 
We  consider  modulation  at  the  10.6  pm  C02  laser  wavelength.  At  this 
wavelength  the  CdGeAs2  crystal  is  transparent  without  any  absorption 
due  to  intraband  transitions  between  the  valence  bands .  With  the 
values  for  the  effective  masses  in  Table  II,  we  find 

{(An)  «  -  AN  X  5-5  X  IQ"19  (cm5)  •  C6-80) 

The  necessary  change  in  the  birefringence  to  change  the  electric  field 
polarization  by  90°  in  a  crystal  with  length  I  is  given  by 

6(An)  -  £  •  (-6M) 

This  corresponds  to  100*  modulation  when  the  crystal  is  inserted  be¬ 
tween  crossed  polarizers.  For  a  2  mm  crystal  the  necessary  injected 
carrier  concentration  is  AN  -  5  x  1015  cm'5  .  This  technique  is  a 
novel  way  to  modulate  a  <X>2  laser.  It  requires  a  material  with 
light  and  anisotropic  effective  masses  such  as  CdGeASg. 
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B.  EXPERIMENTAL  RESULTS 


1 .  Indices  of  Refraction  and  TuninR  Curves 


(a)  CdGeAs,. 

We  have  measured  the  indices  c~  refraction  out  to  10.6  nm  using 
a  L  mm  long  by  3  ram  high  prism  with  an  apex  angle  of  13  degrees.  The 
prism  was  cut  from  boule  No.  16.  This  boule  was  grown  with  0.5  Wt$ 
excess  arsenic.  Table  XIII  lists  the  indices  of  refraction  for  p-type 
material.  CdGeAs.,  has  a  positive  birefringence  of  about  0.1.  For 
the  3.39  urn  and  the  10.6  urn  points  we  used  laser  sources.  The  other 
points  were  taken  with  a  globar  light  source.  Because  of  the  small 
prism  size,  the  amount  of  refracted  light  from  the  globar  was  too 
small  for  the  thermocouple  detector  at  wavelengths  longer  than  five 
microns . 

A  computer  fit  to  the  Sellmeier  expression  in  Eq.  (6.60)  gives 


for  the  ordinary  and  jxtraovdinary  indices.  In  these  expressions  the 


V 


TABLE  XIII 

MEASURED  INDICES  OF  REFRACTION  FOR  CdGeAs2 


n 

e 


n 

o 


long  wavelength  resonance  was  fixed  at  36  ,i®  which  is  the  measured 
reststrahlen  frequency.  Figure  22  shows  a  plot  of  the  indices  of 
refraction.  We  estimate  the  low  frequency  dielectric  constant  from 

Eq.  (6.61)  and  find 


It  is  also  important  to  know  the  temperature  dependence  of  the 
dispersion  and  the  birefringence.  If  the  phasematching  condition  is 
temperature  sensitive,  it  is  possible  to  tune  without  crystal  rotation 
avoiding  annoying  alignment  problems.  For  some  applications  the  tempera 
ture  sensitivity  may  prove  a  disadvantage.  When  the  crystal  has  a 
slight  absorption,  the  resulting  nonuniform  temperature  distribution 
may  make  it  impossible  to  satisfy  the  phasematching  condition  over 
the  whole  crystal  length.  We  have  measured  the  temperature  dependence 
of  the  birefringence  of  CdGeAs,,  at  309  M*  and  10.6  pm.  The  crystal 
was  inserted  into  an  oven  between  crossed  polarizers,  and  the  incident 
light  was  polarized  at  K?  to  the  c-axis.  By  counting  tne  fringes 
in  the  detected  signal  we  determined  the  temperature  dependence  of  the 
birefringence.  With  a  temperature  raise  AT  resulting  in  m  fringes, 

obtain 

d (An)  \m  .  An  a  ,  (6.85) 

dT  i  AT 

where  in  =  n  -  n  Is  the  birefringence,  <  L  the  erystsl  length. 
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and  a  is  the  a-axis  linear  thermal  expansion  coefficient  (cf.  Fig.  10). 
At  10.6  |.un  Of  gives  approximately  a  2$  correction.  The  crystal  we  used 
was  2.96  mm  long  and  it  was  cut  from  boule  No.  28.  Table  XIV  lists 
the  measured  .."'suits  and  Fig.  23  shows  how  the  birefringence 

changes  with  temperature  at  3*39  H1"*  We  have  also  measured  the  tempera- 
ture  dependence  of  the  absolute  indices  of  refraction  at  3.39  4m.  We 
heated  the  prism  slightly  above  room  temperature  And  measured  the  change 
in  the  indices  as  the  prism  cooled  off.  These  results  are  also  listed 
in  Table  XIV.  We  see  from  the  table,  that  the  temperature  dependence 
of  the  birefringence  is  very  large  for  CdGeAs^  and  it  comparable  to 
the  values  for  LiNbO^.^ 

TABLE  XIV 


TEMPERATURE  DEPENDENCE  OF  THE  INDICES  OF  REFRACTION  FOR  CdGeAs2 


X 

(pm) 

T 

(°C) 

dne/dT 

do'1*  V1) 

dnQ/dT 

do'1*  V1) 

d(An)/dT 

do'1* 

Method 

10.6 

-125 

O.389 

fringes 

3-39 

100 

0.772 

fringes 

3-39 

170 

O.833 

fringes 

3-3° 

35 

5-3 

n 

-  1 

prism 

As  discussed  earlier,  the  presence  of  free  carriers  perturbs  the 
indices  of  refraction.  The  prism  we  used  for  the  indices  of  refraction 
measurements  was  p-type,  but  we  did  not  measure  the  carrier  concentration. 
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It  was  probably  between  LO1^  and  LO16  cm  Equation  (6.70)  gives 
the  induced  birefringence  due  to  free  carriers.  With 

2 

e  -15 

r  _  - -  =  2.818  X  10  J  cm 

0  U7Te0mc 


we  find 


^carrier 


(6.86) 


Since  m^  >  m,|  both  for  the  conduction  band  and  the  top  valence  band, 

the  birefringence  decreases  with  increasing  carrier  concentration  both 

.  .-16  -5 

in  n-  and  p-type  material .  For  a  hole  concentration  of  10  cm 
we  have  ,  =  -0.004  at  10.6  pm.  The  indices  are  also  affected 

by  the  Intraband  transitions  between  the  valence  bands,  and  it  is 
important  to  keep  these  effects  in  mind  when  comparing  the  measured 


indices  of  refraction  of  different  CdGeAs2  crystals. 

In  CdGeAs_  phasematched  second  harmonic  generation  is  possible 
between  5  and  18  pm  for  type  I  phase  matching  and  between  5-4  and  15  pm 
for  type  II.  Figure  2b  shows  a  plot  of  the  phasematching  angle  versus 
wavelength.  For  doubling  with  the  SHG  crystal  inside  the  laser  cavity, 
type  I  phasematching  is  most  useful  since,  in  that  case,  the  fundamental 
wave  is  polarized  along  one  of  the  optical  axis.  For  type  II  this  is 
not  possible  and  the  crystal  birefringence  causes  polarization  rotation. 

Since  the  walk-off  angle  limits  the  maximum  useful  interaction 
length,  [cf.  Eq.  (6.15)1,  the  most  attractive  parametric  oscillator 
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construction  would  be  a  90°  phasematched  oscillator  with  a  tunable 
pump.  Figure  25  shows  the  calculated  tuning  curve  for  such  a  device. 

Parametric  tuning  by  crystal  rotation  is  possible  over  almost  the 
whole  transmission  range  of  CdGeAs^.  Figures  26(a)  and  26(b)  show 
tuning  curves  for  several  pump  wavelengths.  A  single  crystal  can  scan 
approximately  20  degrees .  The  walk-off  angle  for  the  extraordinary 
wave  is  typically  about  one  degree.  The  curve  for  Xp  -  4  pm  in 
Fig.  26(a)  represents  an  interesting  special  case  where  parametric 
tuning  is  possible  between  6  and  13  urn  without  changing  the  crystal 
position.  The  large  phasematching  bandwidth  might  be  useful  in  short 
pulse  work. 

A  doubled  C0^  laser  is  a  useful  pump  source  for  CdGeAs^  since  the 
5.3  pm  wavelength  is  slightly  less  than  half  the.  bandgap  frequency  and 
therefore  avoids  two -photorr  absorption.  Figures  27(a)  and  27(b)  show 
the  calculated  tuning  curves  for  the  two  phasematching  conditions. 

The  bandwidth  over  most  of  the  tuning  range  varies  between  10  and  20  cm 
for  a  1  cm  crystal.  For  type  II  phasematching  It  is  possible  to  tune 
through  the  degeneracy  point  and  to  make  a  choice  between  fast 
tuning  and  large  bandwidths  or  slow  tuning  and  small  bandwidths.  Note 
that  the  curvescover  the  very  important  8  to  13  pm  atmospheric  window. 

(b)  CdGePg 

Table  XV  lists  the  measured  indices  of  refraction  for  CdGeP^. 

We  used  a  prism  3.5  mm  by  2  mm  with  an  apex  angle  of  14°.  The  crystal 

5 

was  n-type  and  had  a  resistivity  of  approximately  10  ft-cm.  We 
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WAVELENGTH  (^nt) 


PHASE  MATCHING  ANGLE  (deg) 


FIG.  26a--Theoretical  tuning  curves  for  type  I  phasematching 
in  CdGeHSg  for  several  pump  wavelengths . 
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FIG.  26b _ Theoretical  tuning  curves  for  type  II  phasematching  in  CdGeAs^  for  several 

pump  wavelengths. 
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FIG.  27(a) --Theoretical  tuning  curve  and  minimum  bandwidth 
for  type  I  phaaematching  for  a  one  cm  CdGeAa2 
crystal  pumped  by  a  pump  wavelength  of  5*5  pm. 
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PHASE  MATCHING  ANGLE  (deg) 

FIG.  27(b) — Theoretical  tuning  curve  and  minimum  bandwidth  for  type 

II  phasematching  for  a  one  cm  CdGeASg  crystal  pumped  by  a 
pump  wavelength  of  5.5 
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MINIMUM  BANDWIDTH  (cm1) 


TABLE  XV 


Wavelength 

X((im) 

n 

e 

n 

0 

n  -  n 
e  0 

0.85 

3.U83U 

3.1*283 

0.0551 

0.90 

3-1*397 

3.3958 

0.01*39 

0.95 

3.1*099 

3.3710 

O.0389 

1.00 

3.581*1 

3.3l*93 

O.O3U8 

1.15 

3.3331* 

3.3068 

0.0266 

2.10 

3.2282 

3.211*7 

O.OI35 

3-39 

3.2026 

3.1913 

0.0113 

3.1963 

3.I856 

0.0107 

5.00 

3.1899 

3.1793 

0.0106 

6.00 

3.1853 

3.171*6 

0.0107 

10.6 

3.1517 

3.1367 

0.0150 

estimate  the  index  values  to  be  good  to  at  least  one  part  in  the  third 
pi  ace.  The  average  index  of  3. 2  agrees  well  with  3*1  as  predicted 
from  III-V  analogs ,  but  is  smaller  than  the  value  of  3.6  reported 
by  Goryunova  et  al.^1  Figure  28  shows  a  plot  of  the  indices  of 
refraction.  The  birefringence  of  0.01  is  too  small  for  parametric 
interactions.  The  result  is  in  complete  disagreement  with  Goryunova 
et  al.,  but  agrees  with  Boyd  et  al.^  Goryunova  et  al.  did  their 
measurements  on  platelets  grown  by  vapor  transport  while  we  used 
Bridgman  grown  material.  Assuming  the  Russian  measurements  are  correct, 
the  only  explanation  may  be  that  there  exists  two  forms  of  CdGePg 
both  with  chalcopyrite  structure  but  with  different  order  of  the  Cd 
and  Ge  atoms  within  the  unit  cell.  This  needs  further  clarification 
by  careful  x-ray  analysis. 

The  birefringence  of  CdGePg  similar  to  CdGeASg,  increases  with 

increasing  temperature.  At  2.10  pm  we  measure  dng/dT  =  2.5  X  10  /°C 

and  dn  /dT  =  2.2  X  10"V°C  which  gives  d(An)/dT  «  0-3  X  10  V°c  In 
o' 

the  temperatv.e  range  of  25°C  to  U5°C. 

2 .  Measurement  of  the  Second  Order  Nonlinear  Coefficient 

We  have  measured  the  nonlinear  coefficient  relative  to  GaAs  by 
second  harmonic  generation  of  a  Q-switched  C0,3  laser  using  the  wedge 
technique  described  in  Section  A. 4.  Figure  2 9  shows  a  schematic  of 
the  experimental  set-up.  We  Q-switched  the  C0g  laser  by  rotating  the 
grating  A  He-Ne  laser  reflecting  off  the  grating  provided  a  trigger 
signal.  By  adjusting  the  aperture  the  laser  operated  in  a  Gaussian 
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FIG.  2Q --Measurement  of  the  nonlinear  coefficients 


mode  with  a  pulse  length  of  about  500-400  nsec  and  peak  power  of 

approximately  200  watts.  A  pyroelectric  detector  monitored  the  C0o 

laser  output.  The  detector  was  sufficiently  fast  to  completely  resolve 

the  pulse.  An  uncoated  4  cm  germanium  lens  focused  the  laser  to  a 

40  (jra  spot  size  at  the  wedged  SHG  crystal.  Another  germanium  lens 

focused  the  S!1G  signal  onto  a  liquid  nitrogen  cooled  InSb  detector 

preceded  by  a  sapphire  filter  which  blocked  the  10.6  um  CO,-  laser 

2 

radiation.  The  detector  had  an  active  area  of  0.5  X  0.5  mm  and  a 
response  time  of  8  (isec  so  it  did  not  resolve  the  pulse.  The  detector 
output  was  fed  into  a  boxcar  integrator  and  a  chart  recorder. 

The  CdGeASg  samples  for  the  nonlinear  coefficient  measurements 
came  from  boule  No.  16.  The  GaAs  reference  sample  was  supplied  by 

Q 

Monsanto.  It  was  Cr  doped  and  had  a  resistivity  of  3  x  10  fl-cm. 

Table  XVI  lists  the  parameters  for  the  experiment.  The  measured 
coherence  length  for  CdGeAs^  of  22  *  1  um  agrees  well  with  the  21.2 
calculated  from  Eq.  (6.62)  and  the  index  data.  The  coherence  length 
for  GaAs  is  in  good  agreement  with  the  published  values  in  references  11 
and  62  of  104  ±  7  pa  and  107  -  5  respectively.  For  the  spot  size 

and  crystal  thickness  we  used  both  |  and  B  «  1  [cf.  Eq.  (6.10)]. 

In  addition,  since  the  spot  size  was  much  smaller  than  the  translation 
distance  between  two  SHG  minima,  we  can  use  Eq.  (6.65)  to  determine  the 
nonlinear  coefficient.  The  CdGeAs2  nonlinear  coefficient  measured 
relative  to  GaAs  is 


d  g(CdGeAs2) 


*  3.4  ±  20$ 


(6.87) 


TABLE  XVI 


.-21 


62 


and  with  d^GaAs)  =  1.14  x  10  (mks),  we  have  that  d^g(CdGeAs0)  = 


*14 

3.9  X  10"" 1  (inks) . 


The  optical  properties  of  the  II-IV-V^  compounds  can  be  estimated 

from  the  III-V  analogs.  Using  the  definition  in  J!q.  (2.2)  we  compare 

the  Miller’s  6  for  CdGeAs^  with  the  III-V  analogs  GaAs  and  InAs. 

We  obtain  6(CdGeAs2)/6(GaAs)  =  1.8  and  6(CdGeAs2)/6(lnAs)  =  1.3 

Assuming  that  the  average  bond  nonlinearity  for  CdGeASg  is  similar  to 

l/2 (GaAs  +  InAs)  we  expect  Miller's  6  for  CdGeAs^  to  be  the  average 

of  the  III-V  analogs.  The  measured  value  is  somewhat  larger  but  it 

is  much  smaller  than  calculated  by  Chemla.^  He  uses  Levine’s  bond- 

charge  model  and  obtains  6(CdGeAs2)/fi(GaAs)  =  5.1u 

29 

Recently  Boyd  et  al.,  have  measured  the  nonlinear  coefficient 
for  CdGeAs2.  They  report  d^g(CdGeAs2)/d1^(GaAs)  =2.62  ±  15$  ,  and 
for  the  coherence  length  they  measure  21.5  ±  1  pm  and  calculate  23-0  pm. 
These  results  are  in  good  agreement  with  our  previously  published 
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results.  Also  Goryunova  et  al.,  have  measured  the  nonlinear  coef¬ 
ficient  for  CdGeAs_.  They  used  a  ruby  laser  with  both  the  fundamental 

C 

and  the  second  harmonic  frequency  well  above  the  bandgap  f  requency . 
Because  of  the  reduced  electronic  contribution  to  the  nonlinear  coef¬ 
ficient,  they  measured  a  smaller  value  than  expected  in  the  trans¬ 
parent  region  of  the  cryptal. 

In  Section  VI.A.l(b),  we  defined  three  figures  of  merits  where 

M^  is  applicable  for  90°  phasematching.  For  phasematching  at  other 

angles  M_„_^  ^  or  M_„_  must  be  used.  They  are  the  figure  of 

b  SHGtot  SHGgaus  J  b 


/ 
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merits  for  maximum  SHG  conversion  and  for  SHG  into  a  Gaussian  mode. 

MSHGgaus  also  aPPlies  when  crystals  are  compared  for  use  in 

parametric  oscillators.  Table  XVII  compares  CdGeAs2  with  tellurium 

and  Ag^AsSj  which  are  two  of  the  most  popular  infrared  nonlinear 

materials.  For  reference  we  also  list  GaAs .  CdGeAs2  has  the  second 

largest  known  nonlinear  coefficient  of  any  phasematchable  crystal. 

Only  tellurium  has  a  larger  coefficient,  but  since  tellurium 

also  has  a  larger  double  refraction  angle,  the  C02  laser  doubling 

efficiency  is  largest  in  G'dGeAs2  .  The  maximum  useful  crystal  length 

^max  for  CdGeAs2  is  2.7  mm  for  confocal  focussing  compared  to  only 

0.1  mm  for  tellurium.  The  theoretical  10.6  urn  conversion  efficiency 

in  CdGeAs2  for  a  2-7  mm  crystal  and  confocal  focussing  is  P  /p 

SH'  Fund 

°-001  PFund^  for  the  two  Phasematching  conditions.  This  neglects 

surface  reflections.  We  use  the  calculated  doubling  efficiency  to 

estimate  the  threshold  for  a  degenerate  parametric  oscillator  pumped 

at  5.3  pm.  Assuming  the  equivalent  loss  due  to  the  finite  build-up 

time  is  negligible  compared  to  the  other  losses  in  the  cavity 

[cf.  Eqs.  (6.34)  and  (6.36)],  we  find  from  Eqs.  (6.32)  and  (6.33) 

that  the  threshold  pump  power  P  is 

th 

pth  =  asai  x  l03(w)  (6.88) 

for  a  doubly  resonant  oscillator  (DRO)  and 

Pth  -  “s  X  lo3<W>  (6.89) 


- 135  - 


TABLE  XVII 


FIGURE  OF  MERITS  FOR  INFRARED  NONLINEAR  CRYSTALS 


CdGeASg 

Te 

— 

Ag^AsS^ 

GaAs 

n 

3.60 

4.86 

2.52 

3-3 

d  (rel.  GaAs) 

3.h 

6.9 

o.l4 

1 

Mx  =  d2/n3 

193 

322 

l 

22 

0  (SHG  of  10.6  (am) 
m 

35°31+/(i) 

54°51  (11) 

i4°lo' 

22°30/ 

— 

p 

1°20' 

5°42' 

3°30/ 

— 

'max  * 

2.72 

0.11 

0.56 

— 

"sHGtot  ’  "l 

607 

274 

1 

— 

MSHGgaus  =  Ml/p 

1330 

121 

1 

— . 
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Cor  a  singly  resonant  oscillator  (SRO).  With  10$  signal  and  idler 
losses,  the  threshold  pump  power  and  power  density  are,  respectively, 


Pth 

=  10  W 

and  p^/A  =  1  MW/cm^ 

(6.90) 

for 

a  DRO  and 

Pth 

=  200  W 

and  p  /A  =  20  MW/cm^ 

(6.91) 

for 

a  SRO.  The  area 

A  is  as 

usual,  defined  as  7rw  /2 

P' 

.  Equations 

(6.90)  and  (6.9I)  are  for  confocal  focusing  in  a  2.7  mm  crystal.  A 
larger  crystal  does  not  increase  the  gain,  but  it  reduces  the  power 
density  such  that  crystal  damage  is  less  probable.  The  necessary  pump 
power  to  reach  thiashold  is  available  at  wavelengths  near  5  either 
from  a  TEA  CO  laser  or  from  a  CdGeAs^  doubled  TEA  CO^  laser.  However, 
the  present  CdGeAs^  crystals  are  too  lossy  to  attempt  a  parametric 
oscillator  experiment.  Assuming  the  losses  will  reduce  with  further 
crystal  development,  a  2.7  mm  crystal  should  be  adequate  for  a  doubly 
resonant  oscillator.  For  a  singly  resonant  oscillator,  however,  the 
crystal  should  be  at  least  5  111111  to  reduce  the  power  density.  The 
measured  optical  damage  threshold  for  CdGeAs^  is  40  MW/cm  at  10.6  um. 

3.  Phasematched  Parametric  Interactions 

We  have  demonstrated  several  phasematching  experiments  in  CdGeAs^. 

The  measured  phasematching  angles  are  in  reasonable  agreement  with  the 
calculated  values  based  on  the  measured  indices  of  refraction.  Table  XVIII 
lists  the  results.  For  the  SHG  experiments  we  used  a  CW  CO^  laser. 
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TABLE  XVIII 


PHASEMATCHING  EXPERIMENTS 


Process 

SHG  of 
10.6  pm 

SHG  of 
10.6  pm 

SHG  of 
10.6  pm 

Mixing  of 

10.6  pm  +  5.3  pm 

Boule  No. 

16 

B4 

31 

31 

Phasematching 

condition 

II 

II 

I 

I 

Measured  phase- 

1*9°  ^0' 

48°  50' 

31°  40' 

55°  43’ 

matching  angle 

±  15' 

±  15' 

±  20* 

±  lc 

Calculated  phase¬ 
matching  angle 

5l4°  51’ 

54°  51 ’ 

35°  34’ 

35°  30’ 

Walk-off  angle,  p 

1°  20' 

1°  20' 

1°  20’ 

Crystal  length 
l  (mm) 

1.34 

4.2 

1.95 

1.95 

Measured  angular 
half-width  (deg) 

7- 5-8. 5 

2.25-5.0 

5 

1.7-1 .8 

Calculated  angular 
half-width  (deg) 

8.5 

2.25 

5 

1.9 

We  chopped  the  beam  to  reduce  the  crystal  hewing.  A  motor  scanned 
the  crystal  through  the  phasematching  position.  Figure  30  shows  the 
second  harmonic  output  power  versus  the  external  rotation  angle  for 
one  of  the  crystals.  It  shows  the  regular  sinc-curve  dependence  given 
by  Eq.  (6.20).  The  curve  has  some  oscillations  near  the  maximum. 

These  occur  when  the  laser  is  normal  to  the  crystal.  The  feed-bav.k 
then  makes  the  laser  unstable  which  allows  an  accurate  measurement  of 
the  phasematching  angle.  The  oscillations  are  monitored  at  the  SHG 
wavelength  or  by  recording  the  CO,  laser  output.  The  angle  between 
the  center  of  the  oscillations  and  the  maximum  of  the  sinc-curve 
gives  the  phasematching  angle  relative  to  the  normal  of  the  crystal 
face,  and  by  x-ray  diffraction  we  determine  the  angle  between  the 
normal  and  the  c-axis.  With  this  technique  we  measure  the  phase- 
matching  angle  within  15  minutes. 

In  Table  XVIII  the  external  half-angles  at  the  half-power  points 
are  calculated  using  Eqs.  (6.28),  (6.29),  and  (6.30).  By  comparing  the 
calculated  with  the  measured  half-angles,  we  obtain  information 
about  the  crystal  quality.  In  nonuniform  crystals  with  variations  in 
the  indices  of  refraction,  the  phasematching  condition  may  not  b^  satis¬ 
fied  over  the  whole  crystal  length.  This  leads  to  a  larger  angular 
half-width  than  calculated.  The  crystals  listed  in  Table  XVIII  have 
lengths  between  2  and  U  millimeters  and  the  quality  is  good  enough  for 
the  phasematching  length  to  be  approximately  equal  to  the  crystal  length. 
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i  j  nc  thp  second  harmonic  output  power 

FIG.  30— Angular  dependence  o£  tne  sev-um. 
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To  further  test  the  crystal  quality,  the  crystal  is  translated 
perpendicularly  through  the  laser  beam  while  the  phasematched  second 
harmonic  output  is  monitored.  Figure  Jl  shows  a  typical  horizontal 
and  vertical  scan  for  the  SHG  crystal  from  boule  No.  1 6.  "Good" 
crystal  regions  are  over  distances  of  about  1  mm.  It  is  therefore 
possible  to  avoid  the  areas  with  significant  absorption. 

The  mixing  experiment  in  Table  XVIII  between  10.6  nm  and  5*5 
to  generate  third  harmonic  at  5.55  [m  was  achieved  by  first  doubling 
and  then  mixing  in  two  separate  CdGeAs^  crystals.  The  crystals  are 
placed  right  next  to  each  other.  For  the  doubling  we  used  the  4.2  mm 
type  II  phasematching  crystal  and  the  mixing  was  in  a  1*95  ™  crystal 
cut  for  type  I  phasematching. 

We  have  also  generated  third  harmonic  directly  using  a  TEA  CO^. 
laser.  The  calculated  phasematching  angles  for  tripling  are  49°4 ' 
for  type  I  and  6T°0'  for  type  II  phasematching.  No  solution  exists 
for  type  III  phasematching.  The  1.97  mm  crystal  from  boule  B17  was 
cut  for  type  II  phasematching  which  we  expect  to  give  the  largest 
nonlinear  efficiency.  The  crystal  quality  was  not  as  good  as  for 
some  of  the  other  crystals.  Using  Eqs .  (6.27)  and  (6.49)  we 
calculate  an  angular  half-width  of  3°13'  compared  to  a  measured  half¬ 
width  of  4°15 ’ .  This  corresponds  to  a  phaseraatching  length  of  only 
I.5  mm.  The  crystal  could  be  used  for  both  phasematched  type  II  SHG 

and  type  II  THG  and  the  measured  internal  angle  between  the  two  phase- 

o 

matching  positions  of  11.2  ±  1  was  in  good  agreement  with  the  calcu¬ 
lated  difference  of  12.1°. 
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SHG  OUTPUT  SHG  OUTPUT 


The  measured  phasematching  angle  for  the  SHG  cvvstal  from  boule 
No.  16  was  l+9°30'  as  stated  in  Table  XVIII.  At  one  spot  on  the  crystal, 
however,  we  measured  50°U0 1 .  This  wf»s  probably  due  to  a  crystallite 
at  a  slightly  different  orientation  to  the  rest  of  the  crystal.  To 
determine  how  a  change  in  the  dispersion  D  and  the  birefringence  B 
at  the  fundamental  frequency  affect  the  phasematching  angle  9^  ,  we 

differentiate  Eqs .  (2*5)  and  (2.6)  and  find 


60m  =  2  tan  0m 


fjjD  .$B 
i  D  B  , 


(6.92) 


for  both  phasematching  conditions.  One  of  the  problems  with  absorption 

in  the  crystals  is  that  the  resulting  temperature  increase  changes  the 

phasematching  angle.  We  estimate  the  effect  of  a  temperature  raise  of 

10°C  by  using  =  3-9  X  10'5  V1  from  Table  XIV  and  assuming 

i*lx  10“**  v1  .  For  CdGeAs-  with  a  birefringence  of  B  =  O.O89 
dT  2 

and  a  dispersion  of  D  =  0.029  this  gives  60ffl  =  +  0.6°  for  type  I 

and  60  =  1.2°  for  type  II  phasematching  for  SHG  of  10.6  um.  The 

m 

phasematching  angle  therefore  increases  with  increasing  temperature. 

A  change  in  the  carrier  concentration  also  affects  the  phasematching 
angle.  An  increase  of  the  hole  concentration  of  5  X  10  cm  results 
in  a  decrease  in  the  birefringence  at  10.6  um  of  6B  =  -  0.002  accord¬ 
ing  to  Eq.  (6.86).  From  Eq.  (6  68)  we  calculate  the  change  in  the  dis- 

-U 

persion  of  the  ordinary  index  s  6D  =  0.7  X  10  .  This  increases 

the  phasematching  angle  by  =  0-5°  for  type  I  and  6 0m  =  1°  for 

type  II  phasematching.  Finally  it  is  worth  noting  that  the  indices 
of  refraction  for  p-type  material  are  also  perturbed  by  the  intraband 
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transitions  between  the  split  valence  bands.  It  is  more  complicated, 

however,  to  evaluate  the  magnitude. 

29 

Recently,  Boyd  et  al .  have  measured  the  type  I  and  type  II 
phasematching  angles  for  SHG  of  10.6  um  to  be  35  -  1°  ar»d  51*6  -  0-5° 
respectively.  This  is  in  reasonable  agreement  with  our  results  con¬ 
sidering  they  hnd  crystals  with  much  larger  absorption.  A  bire¬ 
fringence  change  from  boule  to  boule  was  observed  by  Boyd  et  al. 
for  ZnGeP„  where  they  report  a  change  of  0.001  out  of  a  total  bire- 
fringence  of  0.0W3  at  1 .06  (.tm.  Unfortunately  they  did  not  measure  the 
carrier  concentration. 

h.  Maximum  10.6  um  Doubling  Efficiency  and  Optical  Damage  Threshold 

We  have  observed  a  maximum  doubling  efficiency  in  an  uncoated 

CdGeAs  crystal  of  l.U#  for  a  TEA  CO^  laser.  If  we  subtract  the 

reflection  losses  at  the  crystal  faces,  this  gives  an  expected 

doubling  efficiency  of  3#  for  a  coated  crystal.  This  compares  to 

58  59 

the  best  results  reported  for  tellurium.  *  '  The  laser  had  a 
Gaussian  output  beam  with  a  pulse  length  of  200  nsec  and  a  repetition 
rate  of  lU  ppa.  A  27  cm  focal  length  germanium  lens  focused  the  beam 
into  the  SHG  crystal.  The  laser  power  was  measured  by  a  CRL  power 
meter  and  the  second  harmonic  by  an  Eppley  thermopile.  The  crystal 
was  1.95  mm  long  and  it  was  cut  for  type  I  phasematching.  The  measured 
absorption  at  the  laser  and  the  second  harmonic  frequency  was  g  = 

0.8U  cm"1  and  ac  ,  *»  3  cm  1  ,  respectively.  According  to  Eq.  (6.21) 

5 
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this  reduces  the  theoretical  SHG  conversion  efficiency  by  a  factor  of 
0.40.  At  the  measured  1.4$  conversion  efficiency  the  incident  laser 
peak  power  was  12.5  kW.  Including  losses,  this  gives  a  theoretical 
conversion  efficiency  of  2.4$  for  a  square  laser  pulse.  The  pulse 
shape,  however,  was  more  Gaussian  and  this  reduces  the  theoretical 
conversion  by  a  factor  of  approximately  2/j  to  1.6^  which  agrees  well 
with  the  measured  conversion. 

We  observed  a  saturation  in  the  SHG  conversion  efficiency  before 
the  optical  burn  density  threshold  was  reached.  This  might  be  due  to 
either  a  thermal  effect  or  due  to  generated  free  carriers  resulting  from 
optical  absorption  between  an  impurity  level  and  the  valence  or  the  con¬ 
duction  band.  With  reduced  optical  absorption  we  believe  the  satura¬ 
tion  effect  will  vanish.  For  the  measured  crystal  the  saturating  laser 
power  density  inside  the  crystal  was  approximately  15  MW/cm*". 

The  optical  bum  density  measured  by  the  TEA  CO^  laser  was 
58  ±  2  MW/cm?  for  most  crystals.  In  one  crystal,  however,  it  was  as 

high  as  55  MW/cm'.  We  have  also  measured  the  cw  damage  threshold. 

2 

At  10.6  um  it  is  smaller  than  1000  W/cm  . 

5.  Carrier  Absorption 

Figure  J2  shows  the  free  electron  absorption  in  n-type  CdGeASg. 

The  5.5  power  dependence  is  characteristic  of  charged  impurity  scat¬ 
tering  [cf.  Eq.  (6.71)],  and  the  result  agrees  with  reference  66 
which  reports  a  power  dependence  of  J.2  -  5-4. 


Figure  shows  the  lowest  measured  absorption  In  p-type  CdGcAs, 
By  writing  c  =  0.151  eV  -  to  In  Eqs.  (6.75)  and  (6.77)  and  adjusting 
the  hole  concentration  to  2.4  x  1015  cm"3  we  obtain  excellent  agree¬ 
ment  between  the  measured  and  the  calculated  absorption  for  transi¬ 
tions  between  the  valence  bands  Vl  and  v£  .  The  measured  valence 
band  splitting  energy  is  therefore  0.15  ±  0.01  eV  which  Is  not  too  far 
off  our  estimated  0.20  eV .  A  theoretical  value  for  the  splitting 
energy  of  0.17  eV  is  given  In  reference  66.  This  reference  also 
discusses  Intraband  absorption  but  uses  a  formula  that  does  not  apply 
for  the  chalcopyrltes . 

In  conclusion,  in  order  to  obtain  negligible  absorption  In  p-type 
CdGeASg  the  hole  concentration  should  be  smaller  than  1015  cm’5. 
Compensation  or  doping  to  reduce  the  carrier  concentration  Is  therefore 
necessary  to  obtain  high  optical  quality  material. 
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CHAPTER  VII 


APPLICATION  OF  TUNABLE  RADIATION  FOR  AIR  POLLUTION  DETECTION 


Perhaps  the  most  challenging  application  of  tunable  infrared 
radiation  is  for  air  pollution  detection.  In  this  chapter  we  compare 
the  detection  sensitivity  using  an  infrared  parametric  oscillator  to 
other  laser  detection  methods.  The  important  results  have  been 
published  and  are  presented  in  Appendix  D.  Of  the  three  analyzed 
schemes:  Raman  backscattering,  resonance  backscattering,  and  resonance 
absorption;  only  the  last  is  sensitive  enough  to  detect  dispersed 
pollutants. 


j 
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CHAPTER  VIII 


CONCLUSION 


CdGeASg  is  very  promising  for  infrared  parametric  generation 
ard  efficient  SHG  and  mixing.  We  have  measured  the  nonlinear  co¬ 
efficient  and  determined  the  phasematching  conditions.  The  crystal 
is  phasematchable  over  most  of  its  transparency  range  between  2.1* 
and  18  (im  and  it  has  the  highest  figure  of  merit  of  any  phasematchable 
nonlinear  crystal.  CdGeP^  is  not  phasematchable.  It  has  approximately 
the  same  tetragonal  distortion  as  CdGeAs^  but  the  birefringence  is  a 
factor  of  ten  smaller. 

The  growth  of  high  optical  quality  chalcopyrite  crystals  is 
hindered  by  several  problems.  Extensive  crystal  cracking  occurs  during 
growth  and  the  boules  have  nonuniform  optical  transmission.  We  have 
measured  the  linear  thermal  expansion  coefficient  for  CdGeAs^  It  is 
very  anisotropic  and  this  probably  explains  the  large  number  of  cracks. 
We  have  been  able  to  reduce  the  cracks  by  growth  of  single  crystals  in 
a  bismuth  solution.  There  has  been  no  correlation  between  the  non- 
uniform  optical  transmission  and  the  impurity  content.  The  material 
probably  undergoes  significant  self-compensation  during  growth  and  we 
suspect  that  slight  stoichiometric  variations  may  explain  the  non¬ 
uniformity  in  the  optical  transmission. 

To  explain  one  of  the  absorption  mechanisms  in  p-type  material 
we  have  determined  the  bands tructure  near  the  Brillouin  zone  center. 
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The  band structure  results  let  us  calculate  the  absorption  due  to  intra¬ 
band  transitions  between  the  split  valence  bands  and  we  obtain  good 
agreement  with  measured  values  for  CdGeAs^.  Experimentally  we  determine 
the  valence  bend  splitting  energy  to  be  AE  =  0.15  -  0.01  eV. 

V  2 

To  have  negligible  intraband  absorption  the  hole  concentration  must 

15  -3 

be  smaller  than  10  cm 

The  free  carriers  also  affect  the  indices  of  refraction  and  the 

measured  phasematching  angles.  An  increase  in  the  carrier  concentra- 
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tion  of  5  X  10  '  cm  changes  the  phasematching  angle  by  approximately 
1°  for  SHG  of  a  CO^  laser  with  type  II  phasematching. 

Because  of  the  anisotropic  effective  masses  in  CdGeAs^  ,  carrier 
injection  introduces  a  birefringence  change.  We  have  proposed  to 
use  this  effect  to  construct  a  10.6  |im  modulator. 

Since  the  carrier  concentration  is  very  important  for  the  measured 
optical  properties  of  CdGeAs^  we  have  determined  the  temperature  depend¬ 
ence  of  the  carrier  concentration  in  p-type  CdGeAs^.  It  is  constant 
from  77°K  to  200°K  and  then  starts  to  increase.  At  room  temperature  the 
carrier  concentration  in  the  CdGeAs^  grown  by  the  vertical  Bridgman 
method  is  usually  between  5  X  10  ?  to  5  X  10  cm  ^  . 

By  SHG,  THG,  and  mixing  experiments  we  have  verified  several  of 
the  calculated  phasematching  conditions.  We  have  observed  as  high  as 
2$  10.6  nm  doubling  efficiency  in  a  2  mm  CdGeAs^  crystal, and  the 
measured  burn  density  is  40  MW/cm2  for  a  200  nsec  C0^  laser  pulse. 

We  expect  the  quality  to  improve  with  further  crystal  development. 
It  is  a  long  and  tedious  process,  however,  to  determine  the  optimum 
growth  conditions.  To  minimize  the  optical  absorption  it  is  necessary 
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to  reduce  the  carrier  concentration  by  compensation  or  doping.  The 
CdGeASg  has  an  extensive  tuning  range  and  due  to  a  large  figure  of 
merit  only  a  few  millimeter  long  crystal  is  necessary  to  reach  parametric 
oscillator  threshold.  Presently  available  crystals  have  almost  low  enough 
optical  absorption  to  allow  the  construction  of  a  doubly  resonant 
oscillator.  The  calculated  threshold  pump  power  at  5.3  ^m  for  a  2-7  mm 
crystal  is  10W  for  10#  loss  at  the  signal  and  the  idler  wavelength. 

The  most  straight  forward  application  of  an  infrared  parametric 
oscillator  is  in  spectroscopy  with  air  pollution  detection  as  a 
particular  stimulating  possibility. 


APPENDIX  A 


TYPE  I  AND  TYPE  II  PHASEHATCHING  EQUATIONS 


T!  i  extraordinary  index  of  refraction  for  a  wave  propagating  at 

.  8 

an  angle  0  with  respect  to  the  crystal  c-axis  is  given  by 

1 

When  the  birefringence  B  =  ne  "  n0  is  much  smaller  than  nQ  >  we 
can  rewrite  Eq.  (A.l)  as 

n  (0)  n  +  B  sin20  •  (A*2) 

e  o 


2 

cos  0 


sin  0 


(A.l) 


n 


The  phaseraatching  conditions  for  second  harmonic  generation  in  a  positive 
birefringent  crystal  (B  >  0)  is 


(A.5) 


for  type  I  and 


(A.4) 
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/ 


for  type  II  phasematching.  We  calculate  the  phasematching  angle  by 
substituting  Eq.  (A. 2)  into  Eqs .  (A. 3)  and  (A.M-  Defining  the 
dispersion  as  D  =  nf>  -  n“  ,  we  find  for  the  two  phasematching 

conditions  that 

sin  eT  =  fijF  (A-5) 


and 


sin  9 


II 


(A.6) 


where  is  equal  to  the  birefringence  at  the  fundamental  frequency  tu  . 

Similarly  for  a  negative  birefringent  crystal  (B  <  0)  where  the 
phasematching  conditions  for  type  I  and  type  II  phasematching  can  be 
written,  respectively,  as 


f 

i 


For  the  extraordinary  wave,  the  power  flow  and  the  wavevector  are 

not  in  the  same  direction.  The  Poynting  vector  walk-off  angle  is 

8 

given  by 


tan  p 


ne(e)2 

2 


,  (A. 11) 


which  with  the  approximation  |  B  |  «  nQ  reduces  to  the  simpler 

expression 


B 

tan  p  -  —  sin  20 

n 

o 


(A. 12) 


Using  Eqs.  (A. 5)  and  (A. 6)  or  Eqs.  (A. 9)  and  (A.10)  we  can 
determine  quickly  if  a  material  has  sufficient  birefringence  for  phase- 
matched  SHG  . 


appendix  b 


SECOND  HARMONIC  GENERATION  IN  A  HEDGED  PLATELET 

Equation  (6.21)  gives  th*  cur 

8  the  s„c  output  povet  from  a  cry>ul  uich  a 

constant  length  /  t#- 

_  •  a'S"TOS  Ga“aala"  —•  with  .pot  slM8  „ 

w2  wi/  v2  at  the  fundamental  and  thp  e  * 

he  second  harmonic  wave- 

ength.  in  order  to  examine  how  a  wedged  smnnl 

edged  sample  as  in  Fig.  (fi.l) 

-edifies  the  SHG  output  power,  we  recast  Eq.  (6.2l)  in  tems  q£  ^ 
Peak  intensities  and  then  average  the  varying  crystal  length  over  the 

transverse  mode  distribution.  With  p  -  T,  2/o, 

n  ,  2<  Plout  "  ^Ctiw  /3)  and 

2out  ”  we  obtain 


l2 .  issi  4 12  ^ 

2  f”2  +  1  f  (nj  +  l)4 

Here  S(/)  i3  expressed  by 


S(i) 


S  (/)  = 


f"v/” ..'ft  •  -J 


2  Ak  I 


(B.l) 


,  (B.2) 


here  i  -  y  tan  a  according  to  Fig  T 

SHP  nn  t  \  8’  54  '  T°  calculflte  the 

SHG  power,  s(l)  must  be  avera 

8  the  transverse  mode  distri¬ 
bution.  As  before,  we  assume  the  laser  has  a  c.  . 

I([1  ,  .  2,  h8S  a  GaussIan  tode  such  that 

‘  ^  C-  *  ft  .  With  the  crystal  dimension.  In  the 
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APPENDIX  C 


EVALUATION  OF  THE  INTRABAND  ABSORPTION  INTEGRAL 


In  this  appendix  we  evaluate  the  integral  in  Eq.  (6.72).  Denoting 
the  integral  by  S  we  have 


S  « 


ho)  d\  . 


(C.l) 


Let  AE  represent  the  energy  difference  at  k  =  0  between  the  v^  and 
the  v  band.  This  gives 


where  the  definitions  of  Am||  and  Ara^  in  Eq.  (C.2)  should  be  apparent 
and  according  to  Table  II  both  are  positive.  The  v ^  -  v£  band  absorption 
is  most  important  since  this  occurs  in  the  middle  of  the  crystal  trans¬ 
parency  range.  For  this  transition  the  matrix  element  in  Eq.  (6.73)  is 
always  of  the  form 

|M|2  =  Akjj  +  Bk£  .  (C.3) 
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The  matrix  element  for  the  v^  -  v^  transition  also  can  reduce 

to  this  form  for  certain  polarizations.  Using  Eq.  (C.j)  we  have  that 

2  2 

the  integral  in  Eq.  (C.l)  depends  only  on  k j|  and  ki  ,  and 

Sm 

it  is  therefore  useful  to  introduce  polar  coordinates  in  the  (k  ,k  ) 

x'  y' 

plane.  After  the  angular  integration  we  obtain 


S  =  2tt 


/  M 


2 


6 (E -  E^  -  ho)  k|  dk| 


,  (C.fc) 


which  further  reduces  to 


,  2  El/kT 

|  e  '  dk ||  ;  Ex  -  Et  =  ho 


(C.5) 


when  we  make  use  of  Eq.  (C.2)  and  carry  out  the  kj^  integration.  In 

2 

Eq.  (C.5)  ki  is  replaced  by 


L  2 

kj^  =  -  k^  -  —  (A E  -  ho)  Am | 

Am,,  ^ 


,  (c.6) 


We  have  to  consider  separately  the  two  cases  Ae  >  ho  and  Ae  <  ho  , 
and  we  refer  to  the  Integrals  as  respectively  and  »  From  the 
considerations  of  the  Integration  path,  it  follows  that 


lm Am  i  « 


1  J  |H|* 


k  min 


,  (C.7) 


-  l6o  - 


jnd 


imAinl  »  E  /kT 

S2  - -  /  jM|2  e  1  dk(|  ,  (C.8) 

h2  0 


where 

k  >mln  -  2/h2)(AE  -  M&V  •  (C.9) 

With  the  definitions 


<T> 

II 

or 

; 

(C.10) 

and 

Anij 

K  —  -  +  — 

) 

(c.ll) 

“ij  mi,l 

we  can  express  the  Integrals  as 


Sx  -  Sg  [1  -  erf  (  / Ke/kT)] 


+ 


ImAmi  kT 


+  BAm 


1 


K 


(C .12) 
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and 


2nAm]_  kl  /  2rr  Am,  kT 


K 


AAm,|  +  BAmj_  2€Ara|_  B 


K 


kT 


X  exp 


1  ^ 

kT  ral,l 


(C.13) 


In  Eq.  (C.12)  erf(x)  is  the  error  function  defined  by 


erf (x) 


(C.lU) 


In  Chapter  VI,  Section  A. 5  we  use  these  integrals  to  calculate  the 
intraband  absorption  in  CdGeASg  . 
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